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We study the coupling of MSSM fields to heavy modes through cubic superpotential interac- 
tions in F-theory SU(5) GUTs. The couplings are calculated by integrating the overlap of two 
massless and one massive wavefunctions. The overlap integral receives contributions from only 
^ , a small patch around a point of symmetry enhancement thereby allowing the wavefunctions to 

^ I be determined locally on flat space, drastically simplifying the calculation. The cubic coupling 

between two MSSM fields and one of the massive coloured Higgs triplets present in SU(5) GUTs 
is calculated using a local eight-dimensional S0(12) gauge theory. We find that for the most 
natural regions of local parameter space the coupling to the triplet is comparable to or stronger 
than in minimal four-dimensional GUTs thereby, for those regions, reaffirming or strengthening 
constraints from dimension-five proton decay. We also identify possible regions in local parame- 
ter space where the couplings to the lightest generations are substantially suppressed compared 
to minimal four-dimensional GUTs. We further apply our results and techniques to study other 
phenomenologically important operators arising from coupling to heavy modes. In particular 
we calculate within a toy model flavour non-universal soft masses induced by integrating out 
heavy modes which lead to FCNCs. 
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1 Introduction 

Although string theory is primarily motivated as a fundamental unified theory because of its 
ultraviolet behaviour, phenomenological model building within string theory often concerns 
only the infrared spectrum. This is a natural first step given the expected hierarchy between 
the string and electroweak scales. However, heavy modes play a crucial role in our under- 
standing of much of the physics which is relevant to the Standard Model and extensions of it, 
for example by inducing higher dimension operators in the infrared. The fact that studying 
such modes explicitly requires a good understanding of the ultraviolet physics means that this 
is one of the subjects where string phenomenology can play an important role. Heavy modes 
are particularly important in the case of Grand Unified Theories (GUTs). For instance, the 
result of gauge coupling unification at the GUT scale is sensitive to threshold corrections 
from heavy modes, and one of the classic constraints on GUTs comes from dimension-five 
proton decay operators that are induced by integrating out heavy modes. 

Studying detailed properties of these fields, such as their wavefunction profile, is typically 
a difficult prospect because of the complicated geometry associated to realistic models of 
particle physics in string theory. String modes can only be concretely studied in simple 
geometries where a world-sheet description is available. Kaluza-Klein (KK) modes, and the 
closely related Landau-level^, are typically difficult to solve for within some complicated 
Calabi-Yau (CY) geometry. However, in some models, and in particular F-theory (or type 
HB) GUTs, many important operators of the theory are associated to only a small patch 
within the full geometry. The extreme example of this are Yukawa couplings, which are 
associated to just a single point in the geometry. Analogous to the Yukawa couplings there 
are triple couplings between heavy modes and massless modes which can be locally studied 
within a small region around a point. Since locally the complicated global CY geometry is 
decoupled and essentially we can work on flat space, many properties of heavy modes become 
accessible. In this paper we use this local approach to study the coupling of heavy modes to 
massless modes through such a triple coupling operator. This is done by solving for the local 
form of wavefunctions of massive and massless modes and calculating their triple overlap. 

The particular operator that we study, coupling one heavy mode to two massless ones, 
plays a key role in GUTs. One of the general features of SU(5) GUTs is that associated to 
the MSSM Higgs doublets there are coloured triplets which complete a GUT representation. 
These modes have to obtain a mass, leading to the so called doublet-triplet splitting problem. 
Similarly, associated to the Yukawa couplings there are also triple couplings between one 
heavy triplet and two MSSM flelds. Once the heavy triplets are integrated out these couplings 
induce dimension-five baryon and lepton number violating operators that lead to proton 
decay. Thus, understanding such couplings and their flavour structure is of crucial importance 
for placing constraints on GUT models. In minimal field-theory SU(5) GUTs these couplings 
are the same as Yukawa couplings and therefore are exactly known. However, in string theory 
GUTs this is not the case and the couplings can be completely different in nature. This means 
that without knowledge of how the triplets couple to the matter fields and in particular 
whether the coupling to the lightest generations is suppressed in a manner similar to the 

^ These are sometimes referred to as gonions in the intersecting brane literature [1]. 
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Yukawa couplings it is not possible to use dimension-five proton stability to constrain model 
building. The primary aim of this paper is to study the nature of the triplet couplings within 
a realistic string setup thereby performing this crucial step in imposing phenomenological 
constraints on string theory GUTs. 

More generally the paper aims to show that much important physics can be extracted 
by similar calculations of couplings to heavy modes. Indeed, in section [5] we present a toy 
model where such a calculation allows to extract Flavour-Changing-Neutral-Current (FCNC) 
terms which, like proton stability, form one of the important observational constraints on 
ultraviolet physics. We also discuss how our calculations apply to string theoretic realisations 
of the Froggatt-Nielsen mechanism for generating flavour structure this being yet another 
mechanism which relies on higher dimension operators. 

Our focus is on local F-theory GUTs [2H8]. Within this framework a 7-brane carrying an 
SU(5) gauge group wraps a 4-dimensional surface S inside a CY four- fold. Other 7-branes 
intersecting this brane are locally modeled by an enhancement of the gauge symmetry over 
loci in S: along complex curves on S the group enhances by rank 1, to SU(6) or SO(IO), 
while on the points where complex curves intersect it enhances by at least rank 2, to SU(7), 
SO (12) or Eq. We are particularly interested in a point of enhancement to SO (12) as it is 
there that the down-type Yukawa interaction is localised. To describe the physics near such 
point we consider an 8-dimensional gauge theory, which is just super Yang-Mills twisted to 
account for the embedding into the CY four-fold (HEKIl], with S0(12) gauge group broken 
down to SU(5) xU(l) xU(l) by a spatially varying Higgs field. Matter localises onto complex 
curves where the Higgs vev vanishes and at the S0(12) enhancement point three such matter 
curves intersect giving rise to a cubic coupling in the 4-dimensional effective theory. This 
coupling can be calculated directly from dimensional reduction of the 8-dimensional theory by 
integrating the overlap of the internal wavefunctions of localised fields. Yukawa couplings are 
calculated by overlaps of wavefunctions of three massless modes [9j and have been extensively 
studied in [ToHlB] (see also [T7H20] in the context of magnetised D-branes). In this paper we 
calculate the wavefunctions for massive modes around an S0(12) point. Similar calculations 
of massive mode wavefunctions for other models were performed in [21ll22pi6j . Once we obtain 
the wavefunctions for massive modes we can study their overlap with massless wavefunctions, 
thereby probing the cubic coupling discussed above. 

An important property of the dimension-five proton decay operator we are studying is that 
it is a superpotential operator. Since in type IIB string theory and F-theory the superpotential 
does not receive a' corrections, and since integrating out massive string oscillator modes 
induces a' corrections, we do not expect the operator to be induced by exchanging massive 
string oscillator modes. Thus, all the relevant heavy modes which participate in dimension- 
five proton decay are captured within the effective gauge theory described above. 

The calculation of the coupling to massive modes at an S0(12) point is only a part of the 
full calculation required to understand dimension-five proton stability. It is therefore worth 
discussing how the present calculation fits within the full picture of dimension-five proton 
stability in F-theory SU(5) GUTs. We begin by reviewing the constraints on the relevant 
operators. The effective superpotential couplings take the schematic form 

+ YfjT'^ {QiLj + U^Dj) + MTT"^ . (1.1) 
Here Y^ and Y^'^ denote the up- and down-type Yukawa couplings with the MSSM superfields 
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expressed in standard notation. The coloured triplets are denoted by T'" and T'^ and their 
associated triple couplings are and Yf-. In minimal 4-dimensional GUTs we have Y^^'^ = 
y-M|2| rpj^g scale M is related to the mass of the triplets and is expected to be at or below 
the GUT scale, Mqut- Integrating out the heavy triplets leads to dimension-five operators 

W D {QiQjQkLi + UiEjUkDi) . (1.2) 

There are a number of diagrams that lead to proton decay and involve these operators]! At 
TeV scale the diagrams involve a 4-point interaction coming from (ll.2p which has two fermions 
and two scalar superpartners, and a loop factor involving wino or Higgsino exchange to turn 
the scalars into fermions (c.f. figure [1]). This results in nucleon decay to kaons primarily 
(due to the need for a strange quark because of the anti-symmetric colour index). To discuss 
the constraints on the operators ()1.2p let us fix M = Mqut and quote limits on Yj^Yjfi for 
different generation indices El The precise constraints depend on a number of factors such as 



Q , Q 




Figure 1: Schematic structure of Feynman diagrams leading to proton decay through the 
dimension- five effective operators ()1.2p . Similar diagrams exist involving right-handed states. 



the soft masses, the size of the ;U-term and tan /3 (with small tan (3 and large soft masses 
giving generally weaker constraints). Instead of going into the details of these studies, we can 
concentrate on the most relevant aspect for the study in this paper: the difference between 
the Higgs and the coloured triplet couplings. In 4-dimensional field-theory analysis these 
couplings are taken as equal and this leads to an approximate bound Mqut/-^ < 10~^ — 10"'^. 
Thus, if coloured triplet couplings are suppressed with respect to Yukawa couplings by a 
factor larger than the 4-dimensional constraint on Mqut/-^) dimension-five proton stability 
constraints can be satisfied for triplet masses of the order of the GUT scale or above. We can 
estimate the relevant parameters involving the down-type triplet couplings in the minimal 

^This relation and the above superpotential may be slightly modified by more complicated theories where 
the cubic couplings arise after fields in non-trivial GUT representations obtain a vev, as proposed for example 
in [23] to fix the GUT mass relations. 

''There are a large number of papers which study nucleon decay in 4-dimensional supersymmetric GUTs. 
We refer to |24H36) for a subset. Note that a number of these papers were using old experimental results on 
the proton lifetime which has since increased by 2 - 3 orders of magnitude. 

''Note that using the results of ^7] that higher dimension superpotential operators are expected to be 
suppressed by the winding scale, and those of [38H40| showing that the winding scale is also the unification 
scale, implies that Mgut is a quite natural suppression scale. Of course the arguments given are simply 
scaling arguments and should not be taken to hold to significant accuracy. Nevertheless a suppression mass 
scale larger than Mqut seems unlikely given that there are always some heavy modes at or below this scale. 
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4-dimensional field-theory analysis, 



Y'^ Y'^ Y^ Y'^ Y'^ 

-Jh. = 1 -J± 10~1 ^ ^ ^ lQ-2 _Jd ^ _ud ^ Q 3\ 

"V^d ' ' \^d ' \^d \rd ^ ' 

^bb ^bb ^bb ^ bb ^ bb ^ bb 

This simply comes from using the 4-dimensional field-theory equahties, for example 1^'^ = Yfg, 
and the measured quark masses and mixings. We will calculate precisely these ratios for local 
F-theory SU(5) GUTs and compare to the above values to see if there is enough additional 
suppression to avoid proton decay or, alternatively, if there is an enhancement thereby making 
proton decay constraints more severe H 



Calculating the parameters in (jl.3p manifestly requires a theory of flavour. The study 
of flavour structures within F-theory GUTs has been an active research area in the recent 
years [6| [T0HT6lBTVH8] . We make use of the theory of flavour first proposed in [10] and 
subsequently elucidated in [TT] - [T^B2l [T6] The structure is such that all three generations 
are localised on a single matter curve and arise from the degeneracy of massless Landau- 
levels in the presence of flux. This theory of flavour, however, requires ingredients which are 
not present in our setup. More precisely, it was shown in [13l[Illll2] that to generate non- 
vanishing Yukawa couplings for anything other than the heaviest generation requires a non- 
commutative deformation of the theory induced by closed string fluxes or non-perturbative 
effects. We discuss this in more detail in section [5] but for now it is sufficient to state 
that this does not affect the calculation we are performing. Unlike Yukawa couplings, the 
coupling of one massive mode to two massless ones is non-vanishing even in the absence of 
the required non-commutative deformation. Turning on the additional, necessarily small, 
such deformation will only perturb slightly our present calculation thereby maintaining its 
validity. Note that the fact that in the concrete setup we are using the Yukawa couplings 
are rank one in generation space while the triplet couplings are rank three highlights the 
fundamental difference between these couplings. 

The computation that we perform is a necessary one to understand a number of phe- 
nomenological issues. Regarding the particular problem of dimension- five proton stability, it 
is worth discussing some alternative solutions that have been proposed within the context of 
F-theory. One way to avoid inducing proton decay is by having a symmetry which forbids 
it. One such candidate symmetry is a U(l) symmetry, which we label as U(1)pq. Such 
a (massive) symmetry has been studied in detail in F-theory GUT models, see for exam- 
ple |49H6HH3ll45| . Although at the GUT level it was found that many models can exhibit 
such a symmetry, it was shown in a series of papers |53p54ll43ll45P58p 59j that the use of hyper- 
charge flux to break the GUT group and induce doublet-triplet splitting is incompatible with 



^It is important to note that in order to suppress proton decay all the ratios in p.3|l must be suppressed 
since, probing a superpotential coupling, we are working in a weak eigenstate basis. A single large coupling in 
the weak basis can lead to several large couplings in the mass eigenstate basis. Thus, to enhance the rate of 
proton decay it is sufficient that only one of the ratios is larger than in minimal 4-dimensional GUTs. Note 
also that some of the couplings may vanish through other selection rules such as additional symmetries or the 
fact that the colour index in the dimension-five operator must be anti-symmetric and so it cannot involve all 
the same generation. In such cases the ratios involving those operators would not be constrained and it would 
suffice to suppress only the other ratios. 

^There are two key motivations for studying this proposal as opposed to say that of [33] which was based on 
a Froggatt-Nielsen mechanism with additional U(l) symmetries. The first is a practical one: it is not possible 
to study the flavour structure of [43] locally near an SO (12) point. The second is that within the structure 
of [33] the relations (|1.3p are always at least as strong as in minimal field-theory GUTs since the suppression 
by the U(l) symmetries acts on the triplets in the same way as on the Ifiggs doublets. 
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such a symmetrylj The precise statement is that the presence of a U(1)pq symmetry neces- 
sarily imphes the presence of exotic non-MSSM states in the massless spectrum. The mass 
of the exotic states is set by the scale at which the U(1)pq symmetry is broken and therefore 
the constraints coming from dimension-five proton decay translate to constraints on the mass 
of the exotic states. The phenomenology associated to different masses for the exotic states 
was studied in detail for a number of models in |45P591I60] . Since the exotic states do not form 
complete GUT multiplets, the most immediate constraints on their masses come from gauge 
coupling unification. The tension between a large exotics mass to maintain gauge coupling 
unification and a small mass to preserve an approximate U(1)pq symmetry implies that it is 
difficult to practically realise the full suppression necessary for dimension-five proton stability 
using such a symmetry alone. The suppression due to a U(1)pq symmetry is additive to that 
studied in this work and therefore whether we find additional suppression or alternatively an 
enhancement of coupling to massive modes can allow for or rule out a number of proposed 
models. 

An alternative possibility for suppressing dimension-five proton decay even without a 
U(l)pQ is keeping the matter curves associated to the up and down Higgs fields in the same 
homology class but still geometrically separated. The interaction between the up and down 
triplets may then be suppressed by their small wavefunction overlap, although explicitly 
studying this would require a calculation of massive wavefunctions similar to that presented 
in this paper. Apart from the fact that this rather complicated setup has yet to be realised 
explicitly, there are a number of phenomenological problems with such a setup. The first is 
that the use of hypercharge flux for doublet-triplet splitting is difficult since it acts in the same 
way on both the Higgs curves. Another problem is that the theory of flavour of |10) is based 
on local geometric symmetries which means that in order to correlate the up- and down-type 
Yukawas, as is required by a realistic CKM matrix, the geometric separation between them 
should be small. Indeed this is one of the primary motivations presented in [52J for a proposed 
point of Eg uniflcation. 

The outline of this paper is as follows. In section [2] we introduce the effective theory that 
we will be using. Following this, in sections [3] and [4] we present the actual calculations of the 
relevant wavefunctions and their overlaps. In section [5] the results and their phenomenological 
applications are discussed in detail. We summarise our findings in section [6l In appendix 
[Al we present the wavefunctions and overlaps for a more general set of background fluxes, 
consisting also of oblique fluxes. In appendix |B] we study in more detail the normalisation of 
the wavefunctions. 



2 The effective theory 

We consider F-theory on an elliptically fibered Calabi-Yau 4-fold X, with the degrees of 
freedom of an SU(5) GUT localised on a codimension-2 singularity. In the infrared those are 
described by a twisted 8-dimensional M = 1 gauge theory, with gauge group G and support 
on M^'^ X S*, where 5 is a 4-dimensional Kahler sub-manifold of X [ll[5]|j In this section 
we describe the 8-dimensional effective theory and how its massless and massive localised 

^As pointed out in [62], the same problem arises for Wilson-line GUT breaking. 

*If 5* is shrinkable (more formally it has an ample normal bundle), the resulting 4-dimensional gauge 
coupling QfGUT can be tuned independently of Mpianck- However we do not necessarily assume this property 
in our analysis. 
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spectrum is calculated. Related computations to the ones that we describe here have also 
recently appeared in [16]. 



2.1 The 8-dimensional effective theory 

For convenience, we arrange the 8-dimensional fields in adjoint valued, 5-valued, 4-dimensional 
M = 1 multiplets 

-A-m = { A ffi, ip jfi, Qffi) , (2-1) 
V = {v,A^,V) . 

The subindices on the fields denote their local differential structure on S. Thus, for instance 
Affi G Clg (8) ad(i-*) where 17^ denotes the space of holomorphic p-forms on S and P is the 
principal bundle (in the adjoint representation) associated to the gauge group G. Here A 
and $ are chiral multiplets with respective F-terms Q and T-L. V is a vector multiplet with 
D-term T>. and (pmn are complex scalars while tpm, Xmn, and rj are fermions. 

The action for the effective theory was given in [5j. For the bosonic components of the 
multiplets it reads lf| 



58d = M: I (fx Tr 



1 , 1 

2 4 

+nAn- au- a-h-g aOa^-G aBa^p 



+2 (^u; A F(i-i) + c^]) P - 2ia; A F^'f A 4°'^^^ + 



(2.2) 



where uj is the Kahler form of S. Our conventions are such that uj is dimensionless, A and 
if have dimensions of mass and the auxiliary fields V, Q and % have dimensions of (mass)^. 
M^< denotes the UV cutoff of the theory. In the weakly coupled type IIB limit this is related 
to the string scale as 

Mt = (27r)-5a'-2 . (2.3) 

Above this scale, corrections to eq. (|2.2p in the form of higher derivative couplings become 
important and keeping only the leading term in the M^, expansion, eq. (|2.2p . is not a valid 
approximation. Thus, in what follows we shall stick to the regime where, at every point of 
X 5, 

{dA), (dip) < Ml (2.4) 



Ideally we would like to dimensionally reduce the 8-dimensional effective action (j2.2p 
on S in order to obtain the spectrum of 4-dimensional fields with masses smaller than M^,. 
However, such a program would require the precise knowledge of the geometry of S, which 
in general is only available for few highly symmetric spaces such as T^, x or P^ (see for 
instance [18]). Alternatively, we can solve the equations of motion in a local patch around 
a particular point of S where the energy density of a set of charged modes localises. This 
approach has been extensively used in recent phenomenological studies of Yukawa couplings 
in F-theory GUTs (see e.g. |10H16 j). In what follows we describe it in detail. 



^This action can be shown to be equivalent to 8-dimensional super Yang-Mills theory with a non-trivial 
Higgs bundle [H]. 
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2.2 Equations of motion for localised fields 



Let us first consider 4-dimensional massless fields. Setting tlie 4-dimensional variations of 
tlie fields to zero, the equations of motion that follow for their internal wavefunctions are [5] 



n - f(2'°) = , (2.5) 

i [if, (p] + 2uj ^ + = , (2.6) 

2iuj ^g - Ba^ = , (2.7) 

-dH + 2u}AdV + gAif-x^i^- i2V2uj A A = , (2.8) 

ujAdAi^ + ^[^,x]=0, (2.9) 

BaX - 2i\/2a; A Sat] -[ip,ip]=0 , (2.10) 

BAi^-V2[^,r]]=0 , (2.11) 

-V2[fi,x]-BAg-^[i^,^P]=0. (2.12) 



where we have also included the equations of motion for the fermionic fluctuations. Eqs. (j2.5p 
and (12.6P are usually dubbed as the F-term and the D-term conditions for the flux, respec- 
tively. 

Generically, the gauge group G is broken by ((/?) and (A) to a smaller subgroup G' C G. 
We are interested in vacua where {ip) and {A) take values in a subgroup Gh = U{1) x . . . x ^7(1) 
belonging to the Cartan of G, 

G^G'xGh (2.13) 

with G' the commutant of G and Gh- In that case eq. (12. 7p can be simply satisfied by 
requiring {(p) to be holomorphic 

Ba{v) = B{ip) + [A, {ip)] = B{p) = . (2.14) 

Note also that to preserve 4-dimensional Poincare invariance we must impose 

(X) = W = = , (2.15) 

which, making use of eqs. (|2.5p - (|2.7p . imply that eqs. (j2.8p and (j2.12p are automatically 
satisfied. 

Modes charged under Gh arrange into Landau levels and are localised around points in 
S where {ip) = {A) = 0. In order to have a description of these modes it therefore suffices to 
consider a local patch around a localisation point. We can take the Kahler form to be given 
in the local patch by the expansion 

a; = - [dzi A dzi + dz2 A d^) + . . . (2.16) 

where the dots denote higher order terms in the two local complex coordinates zi and Z2 ■ Our 
conventions are such that coordinates are dimensionless, dzi denote local vielbein 1-forms and 
the origin of coordinates is at the localisation point. 

Similarly, we expancj^ 

{A) = —^ lm{M^jZidzj)Qa + . . . , (2.17) 

(99) = A'UR^ m^Zi Qadzi A dza + • • • , (2.18) 
^°Note that most generally we could have also considered constant terms in these local expansions. 



8 



where Mf- and are arbitrary numbers related to the quanta of gauge and Higgs fluxes 
and Qa denote the difl'erent Abehan generators of Gh- We have also introduced the standard 
lengths and R±_ in the local patch and its transverse space respectively, measured in 
units, which for simplicity we have taken to be the same in all directions. Note that 
we have chosen to parameterise the dimensionful part of ip with the transverse ("winding") 
scale -Rj_M*. This is not an arbitrary choice but follows from embedding the 8-dimensional 
theory into a 10-dimensional theory where the Higgs would correspond to deformations of 
the 7-brane into the normal directions. In the 8-dimensional theory the Higgs kinetic term 
arises from the pull-back of the 10-dimensional metric normal to the brane. 



(2.19) 

where ^33 is the metric transverse to the 7-brane and (j)^ is the complex scalar parameterizing 
geometric deformations of the 7-brane along the holomorphic normal vector. In order to have 
a canonically normalised quasi-topological term (c.f. eq. (j2.2p l. in the last equality we have 
redefined, 

if = M^R±L^Q = M^R^(t)^dz^ A dz^ , (2.20) 

with Vt the local holomorphic 3-form of the 3-fold base -B3 of X. Hence, the factor R±M^: is 
the appropriate one for the canonically normalised Higgs field in the 8-dimensional theory^ 



Plugging the above local expansions into eqs. (I2.9p - (l2.1ip we obtain that the relevant 
F-term equations for a given massless 4-dimensional fermionic field read, to leading order in 
the coordinates, 

ID* = , 



with. 



/ 

-D, 
-D2 
\-D3 





-Dl 



D 



D2 
Dl 



-D 



D3 \ 



-D 



and 



D, = 



R 



D3 = -M^RjL Qam'^Zi 



D] 



/ 

' R\\ 



\x J 



(2.21) 



(2.22) 



1,2 



Dl = M,R^qa{m1 



Zi 



(2.23) 
(2.24) 



In these expressions q is the vector of G/^-charges for the localised mode and we have relabeled 
LfYi — )• i-p and X12 — ^ X to simplify the notation. To obtain a finite set of solutions, we have to 



^^The above scalings with K\\ and Kx_ can also be understood from the T-dual setup with magnetised D9- 
branes. Indeed, T-duahsing along the transverse space to the 7-brane, the Higgs and gauge fluxes are mapped 
respectively to gauge fluxes Fia and Fij on a stack of D9-branes [TS]. In a vielbein basis the components of 
the flux are respectively, 

mMl „ MMl 



Fr, 



T-dualising along the transverse directions, these become the Higgs and gauge fluxes on the 7-brane 

MMl 



where R^^ = {R-T) \ in agreement with (PTT)! and (HHH}. 



9 



supplement these equations with a set of boundary conditions encoding the global obstruction 
from the topology of S and which, in particular, determine the degeneracy of the zero modes. 

Massive 4-dimensional fields can similarly be accounted for by the 8-dimensional effective 
theory. In that case, one obtains the more general set of equations 

©"^B^ = ImApl' , (2.25) 

where mx is the mass of the 4-dimensional field and the same definitions above hold. As 
we explicitly show in next subsection, these are the equations of motion for a set of three 
complex quantum harmonic oscillators which can be solved by means of standard techniques 
in quantum mechanics. 



2.3 Localised fields and supersymmetric quantum mechanics 



Let us first solve the equations of motion for 4-dimensional massless fields, eq. (|2.2ip or 
equivalently ec|_. (I2.25P with mx = 0. For that we closely follow the techniques developed 



m 
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Different matter representations have different charges under the gauge group generators 
and therefore different equations governing their wavefunctions. We take ^ to transform in a 
representation TZ of the gauge group, with Di the corresponding gauge covariant derivatives 
defined in eq. (I2.23p . From eq. (|2.22p we observe that the operator which appears in the 
left-hand-side of eq. (|2.25p can be written as 

D'^D = -AI + B , 



where 



and 



A = 



E 

=1,2,3 



D]Di 



(2.26) 
(2.27) 



/O 






[Di, 
[Di, 



dI] [D2,Dl] [DlD2] + [DlD^]j 



it 



it 



(2.28) 



We have made use of the F-term equations for the background, that we take to preserve N >1 
supersymmetry in 4-dimensions, in order to simplify the hermitian matrix B. Imposing also 
the D-term condition on the background it is easy to check that B is traceless. 

A suitable approach to obtain the zero mode wavefunctions is therefore to make a change 
of basis which diagonalises B and to solve the equations of motion in that basis. Let J be the 
matrix which diagonalises B and has canonically normalised column vectors, 

2 



Ri 



— diag(0,Ai, A2,A3 



(2.29) 



where Ai + A2 + A3 = 0. The dimensionless eigenvalues Ap are given by the three roots of the 
characteristic polynomial of the non-trivial part of (-R||/M=k)^B, which is a depressed cubic 
equation. We can rotate the operator D to the diagonal basis by taking. 



•JJ . 



(2.30) 



^^See also [B3] for related work. 
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Notice that B has again the same structure as in (j2.22p but in the new basis covariant 
derivatives are given by 

3 1 ^ 

bp = ^ IkpDk = jT— ^ Cp, kDk , (2-31) 

k=i ll^pllfc=i 

with the p-th eigenvector of B and ||^p|| its norm. In particular, it is simple to check that 
the only non-vanishing commutators of the rotated covariant derivatives are the diagonal 
ones, 

[bl,bp] = -i^-j^j Xp, p = 1,2,3. (2.32) 

As we have advanced, this is the algebra for the ladder operators of a set of three quantum 
harmonic oscillators. 

Generically we can distinguish four towers of solutions to eq. (12.251) . one per eigenvector 
of B. These four towers can be identified with the four complex fermions of a (broken) 
= 4 supermultiplet. In particular, if the three non-trivial eigenvalues Xp are different from 
zero, there is a massless chiral M = 1 supermultiplet in the 4-dimensional spectrum, which 
corresponds to the ground state of one of the above four towers. To see this more explicitly, 
let us assume for a while that A1A2A3 < with two positive and one negative real eigenvalues. 
We take Ai to be the negative eigenvalue. The localised normalisable solution satisfies the 
equations 



We can identify the raising and lowering a operators as, 

ai = ibi , 0,2 = ib\ , 03 = ib\ , (2.34) 



^ 



Diif = 

A 

[dI^ =0 



blip =0 . (2.33) 



a\ = ib\ , a\ = ib2 , a\ = ib^ , 

so that the function ip in eq. ()2.33p is annihilated by the three operators. More generically, 
for fields transforming in the representation IZ, we have in the diagonal basis 

n : ]D)tB= ^ a\ail+[-^\ diag(-Ai, 0, A2 - Ai, A3 - Ai) . (2.35) 

Since the ground state in each of the four towers of fermions is by definition annihilated by 
all lowering operators en, the four entries in the last term correspond to the masses of these 
ground states. Their wavefunctions are given in terms of the function (p as, 

^>p = ^^(zi,Z2,zi,Z2) , p = 0,l,2,3 (2.36) 

where is a normalisation constant. Similarly, wavefunctions for the heavier modes in each 
tower are obtained by acting on the corresponding ground state wavefunction with the raising 
operators. We can label these fields by three quantum numbers, n, m and /, according to 

(ai/M,)"+'+™ - - , 

= , - L12 )"(^2)"'(I^3)'^p , (2.37) 
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where the particular pre-factor ensures the correct normaUsation. Thus, calculating massive 
wavefunctions is a simple task of applying differential operators to functions. The mass of 
the resulting 4-dimensional fields is given by 

Mi,__, = ( ^ ) [-{n + l)Ai + m\, + /A3] , (2.38) 
2 



fO,(n,m,i) \^ ^1 



f2 ^ 



In particular, ^'i, (0,0,0) denotes the wavefunction for the massless chiral fermion transforming 
in the TZ representation. 

Whereas this description is complete for massless chiral fields, massive fields contain 
both chiralities and the above wavefunctions only represent half of their degrees of freedom, 
namely those transforming in the TZ representation of the gauge group. Wavefunctions for 
the TZ components of the massive fields can be worked out following the same procedure, 
taking care of the change of sign in the charges. One may easily check that the analogous 
operator to ()2.35p for fields transforming in the TZ representation is 

1Z : ]D)tD= ^ a\ail+{^\ diag(-Ai, -2Ai, A3, A2) . (2.39) 

j=l,2,3 ^ II ^ 

Wavefunctions are therefore given again by the same functions [nmi)^ with p = 0, 1, 2, 3, 
but the corresponding masses are shifted with respect to eq. (I2.38p . Massive components 
transforming in the IZ and the TZ representations pair up non-trivially. For instance, the first 
excited states of the massless mode in the TZ representation pair up with same-mass ground 
states in the TZ representation^ Note also that there is no massless fermion transforming in 
the TZ representation, as expected. 

Wavefunctions for the scalar fields can be worked out in a similar way and, in particular, 
they can be shown to be identical to the ones of their corresponding fermionic superpartners, 
as a consequence of supersymmetry and fiatness of the local patch. 

Summarizing, we have shown that at each localisation point in S there are four towers of 
fields with equal gauge charges, corresponding to the degrees of freedom of a broken A/" = 4 
supermultiplet. In these conventions, for A1A2A3 < there is a localised massless N = 1 chiral 
supermultiplet transforming in the TZ representation of the gauge group. The degeneracy of 
this field is only globally determined. It is also easy to check that for Ai A2A3 > the roles of TZ 
and IZ are exchanged and there is instead a massless N = 1 chiral supermultiplet transforming 
in the TZ representation. These two possibilities are separated by a wall of marginal stability 
at A1A2A3 = 0. At this wall at least one of the three eigenvalues Ap vanishes and 4-dimensional 
fields arrange into N = 2 supermultiplets (or = 4 supermultiplets if all eigenvalues are 



^■^This non-trivial pairing has its origin in the fact that the 4-dimensionaI mass term comes from the 8- 
dimensional kinetic term, and the latter contains a V^Di operator acting non-trivially on ^ . 
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zero) with conserved Kaluza-Klein momentum. In that case, the wavefunction of the fields 
is no-longer localised along the matter curve, and their mass is determined by the particular 
topology of the curve. We present in section [3] an example of this type. 



2.4 Validity of the local approach 

The wavefunctions and consequently, through their overlaps, the cubic couplings in the 4- 
dimensional theory depend on the parameters of the 8-dimensional theory such as the fluxes 
and the local scales and R±. It is therefore important to quantify the possible range of 
these parameters which is consistent with the local effective theory being used. 

The first constraint we must impose is ()2.4p which ensures that higher derivative correc- 
tions to the 8-dimensional effective action are negligible. Using the expressions ()2.17p and 
(fZTHD this gives 

« 1 , « 1 . (2.40) 

These amount to small intersection angles and small flux densities. The flux parameters M?- 
and would be integer quantised in a homogenous setup but in the local setup need not 
be so. However, generically they are expected to be of order one and we shall therefore take 
them as so while keeping in mind that the local freedom to adjust the fluxes allows for some 
flexibility in satisfying the consistency constraints. Taking the fluxes as such we can rephrase 
(|2.40p in terms of geometric constraints. We define 

R = RnR±, £=^, (2.41) 



^11 



using which we can write 



e < 1 , ^ < 1 • (2.42) 
R 

We can also consider these as constraints on mass scales: it is simple to check that for large 
R the eigenvalues Xp scale as, 

Ai, A2~i?, Ag-l, (2.43) 

and are independent of £. Therefore, from eqs. (12.38P we observe two types of massive modes 
with masses scaling as M ~ M^/R\\ or M ~ M^y/R/R\\. The mass of these modes should be 
kept below the cutoff scale of the theory which is consistent with the constraints (j2.42p . For 
generic order one fluxes the two above constraints can be simultaneously satisfied by taking 

R\\ > , i^ii > 1 . (2.44) 

In this limit there is a large number of 4-dimensional massive fields below the cutoff scale M^, . 
Note that, although not necessarily required, these constraints also allow for length scales 
R± < 1 . The stability of the 8-dimensional effective theory against 1 /Mpianck corrections for 
such small values of R± depends on the particular connection between the local and global 
scales, which we now discuss. 

The relation between the local scales and R± and the global ones is model dependent 
and generically too complicated to be computed explicitly in given models. Whereas the limit 
(I2.44P can always be taken in the local setup, once we begin to relate local scales to global 
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ones new phenomenological constraints are expected to arise from the observed values of 
«GUT and Mpianck- For instance, if S is completely homogenous then we have approximately 




a, 



GUT ■ 



4 



(2.45) 



Since we observe a] 



-1 

'gut 



24 this implies that R\i cannot be too large. The relation is 



approximate and the space S is in general not homogenous, but nonetheless it is difficult to 
conceive a departure of the local scale too far from ()2.45p . We must therefore keep in mind 
that although formally our calculations can be made very precise by taking the limit (I2.44p . 
in a phenomenologically viable setup there will be corrections (essentially a' corrections) 
that are not hugely suppressed. Taking this into account, and that the constraints are only 
approximate up to order one factors, in what follows we allow ourselves to take R in the range 
1 < i? < 25, with e also in the range 1 < < 25 but chosen appropriately such that for 
each value of R eqs. (j2.42p are satisfied. Note that the most natural values are towards the 
lower end of the range, however, due to the strong model dependence of the relation between 
local and global scales this range is only an approximate one and some flexibility should be 
allowed. 

Similarly, we expect the local scale R± to be related to the global ones and in particular 
to Mpianck- The particular relation strongly depends on the geometry of the CY base In 
the case of a torus = we have 



It is therefore important to note that in a torus, and more generally in a near homogeneous 
setup, the observed values of Mpianck and ogut are not compatible with the constraints (|2.42p 
and we expect higher derivative corrections to the effective 8-dimensional theory coming from 
large brane intersection angles. At a deeper level this can be taken as motivation for local 
models based on contractible cycles as then the scaling with respect to the Planck scale 
is expected to be modified to the schematic form R± ~ (Mpianck/-^*)^'^^- More generally 
the geometry can lead to differences between R± and the global scales either coming from 
inhomogeneities of the divisors or from the geometry allowing a decoupling of the intersecting 
brane setup from the overall volume. Given this in general we do not attempt to relate the 
local scales with the global ones. 

3 The SO (12) enhancement point 
3.1 The SO(12) point and background fluxes 

We now apply the procedure described in the previous section to the point in S where the 
down-type Yukawa coupling localises. At that locus there is an S0(12) enhancement of the 
gauge symmetry which can be seen by decomposing the adjoint 



R± 



gs^Planck«GUT 



(2.46) 




) 



(3.1) 
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The spontaneous breaking of the S0(12) symmetry away from the enhancement point can 
be obtained by turning on a background for the Higgs scalar, 



{^) = M,R^ (^Q^ + ^qA , (3.2) 

where vi and V2 are dimensionless parameters. The generators Qi and Q2 are those cor- 
responding to the U(l) factors in the decomposition of S0(12). This Higgs background 
describes three sets of intersecting 7-branes with locaHsed matter on their intersection curves 



10a/ 

5 a/ 
5// 



Z2 = 0, 

zi = 0, 

V2Z1 + ^122 = . (3.3) 



We will be solving for the wavefunctions of modes localised along these curves. 

In order to get chiral matter, U(l) flux must be turned on along the generators Qi and 
Q2. Moreover, to break the GUT group and induce doublet-triplet splitting also flux must 
be turned on along the hypercharge direction of SU(5) so that, 

SU(5) ^ SU(3) X SU(2)^ X U(l)y , 

5 ^ (3,1)1/3 e (i,2)_i/2 , 

10 ^ (3,l)_2/3 ® (3,2)1/6 ® (1,1)1 . 

The geometric properties of these fluxes are dependent not only on S but also on the full 
CY four-fold X. Let us recall some of the defining global properties of the fluxes. First, 
the hypercharge flux must be turned on along a cycle which is homologically non-trivial 
when pulled back to the GUT divisor S, but trivial in the full CY so that U(l)y is massless 
[M1I11E1I3 • Secondly, the fluxes must be turned on such that they induce the correct chiral 
matter spectrum, which is determined as 



«(i,2),/, - ri(i,2)_,/, = M, + N, (3.4) 



for the 5 curves and 



n(3,l)-2/3 - ^(3,1)2/3 = ^^10 - ^ ' 

"-(1,1)1 - "-(1,1)-! = ^^10 + N , (3.5) 

for the 10 curves, where n-ji denotes the number of massless 4-dimensional fields transforming 
in the TZ representation of SU(3)xSU(2)ixU(l)y. Here the fluxes M5, Miq and are 

specified by fractional line-bundles Ly , Vio and V5 such that Miq = deg (^Ly^ (g) Vio^ , M5 = 

deg (^Ly^^^ fgi Vs^ and = deg ^Ly . To maintain complete representations on the matter 
curves we want = for those curves, while to induce doublet-triplet splitting we want 
M5 = and = 1 for the Higgs curve. 

Working locally near a point on S we are not sensitive to the full global structure of the 
fluxes. Indeed, all the geometry of the cycles locally reduces to the four possible components 
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of the flux along dziAdzi, dz2Adz2, dzif\dz2 and dz2f\dzi. These components are constrained 
by the local D-term condition 

^^^(1,1)^0, (3.6) 

whereas the F-term condition simply requires = i^C'^) = Q. 

We assume that locally the flux takes a constant proflle in 5, neglecting a possible spatial 
dependence of the flux. This can be expected to be a decent approximation if the curvature 
around the enhancement point is small, as in that case it can be thought as the leading term 
of a Taylor expansion in the local coordinates, as we have argued in the previous section. 
Taking varying flux into account would lead to technical difficulties, in particular D-terms 
would generally not be solved by a flat-space proflle, which would imply having to work in 
curved space [.14J. Moreover, for simplicity, in the main part of the paper we do not turn 
on flux along the oblique components dzi A dz2 and dz2 A dzi. Oblique fluxes turn out to 
not affect the physics of the wavefunction in a qualitatively important way. We relegate 
a treatment of the more general flux including such components to appendix [A) A general 
suitable choice for the U(l) flux is therefore 



i7(l>l) 



2iMi 



{dzi A dzi - dz2 A dz2){-MiQi + M2Q2 + iQy) 



(3.7) 



where Mi, M2 and 7 are dimensionless real constants and the generator Qy is along the 
hypercharge direction in SU(5)gut- The gauge potential associated with this flux reads 



A 



{zidzi - zidzi - Z2dz2 + Z2dz2){-MiQi + M2Q2 + iQy) 



(3.8) 



The relation between the local values of the flux Mi , M2 and 7 and the global integrated 
values in (|3.4p and (j3.5p is subtle. The connection is that the local fluxes Mi, M2 and 7 
determine the chirality of the localised fields, which generically arrange in M = 1 supermul- 
tiplets as we have described in previous section. At a given localisation point, eq. (j2.2ip 
has an infinite number of solutions. Consistency with the topological data of the fiux and S 
however selects a finite subset of size given by the global integrated values in (|3.4p and (|3.5p , 
in accordance with standard index theorems. 



Using the fiux (13. 7p the chirality for a given localised mode is determined by the analogous 
local expressions to eqs. ()3.4p and (|3.5p . That is, for the 5m curve. 



'^(3,l)-l/3 

'^(1,2)1/2 



sign 
sign 



-Ml + ^7 



-Ml + -7 --7 



(3.9) 



for the 10m curve. 



^(3,2)1/6 

'^(3,l)-2/3 
'^(1,1)1 



Sign 
sign 
sign 



M2 + -7 
6 



M, + -7j--7 
M2 + -7 j + g7 



(3.10) 
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and for the 5h curve, 



%,i)_i/3 = sign 
'^(1,2)1/2 = sign 

where d-ji = +1 {S-ji = —1) means that the corresponding set of locahsed massless 4- 
dimensional fields transforms in the TZ ilZ) representation of SU(3)x 811(2)/, xU(l)y and we 
have taken V2 > vi without loss of generality^ In this regard, it is also worth stressing that 
if the local flux vanishes along the matter curve for a given representation the wavefunction 
does not localise, as we have already commented in the previous section. 

With these relations we see that there are some constraints on the local fluxes that one 
has to satisfy in order to properly model the massless spectrum. In particular, we shall 
require that the expressions in (13. 9p are all negative and that the expressions in (IS.lOp are all 
positive. This can be implemented by taking Mi and M2 positive and much larger than 7. 
On the Higgs curve we require that the second expression of (I3.1ip is negative. The sign of 
the first expression of (jS.lip determines whether locally there is a massless triplet, a massless 
anti-triplet or a vector-like pair, with the mass of the latter depending on the particular 
topology of the matter curve, 

, Massless (3, l)_i/3 (3.12) 
0, Massless (3, 1)1/3 

, Massless or massive (3, 1)1/3 © (3, l)-i/3 • 

We will study all three possibilities in the next sections. The most appealing case is the third 
one, as the mass of such a vector-like mode is only determined globally. If the fiux vanishes not 
only locally but everywhere along the matter curve the vector-like pair is massless whenever 
ho ^S, K^^^^ 7^ 0, where E denotes the matter curve and K its canonical bundle [HIS]. In the 

interesting case where the 5h curve is a on a del-Pezzo such a vector-like pair is massive 
because of the twisting. The other two cases in (j3.12p are different from the third in that 
the only way to give them a mass would be to somehow deform the local geometry. Such a 
deformation could be due, for instance, to the presence of another nearby enhancement point 
with a localised massless triplet of opposite chirality, such that the massless mode (and their 
Landau replicas) develop a supersymmetric //-term. If the shift in the masses is much smaller 
than M*/i?|j we may think of the first two cases in eq. (I3.12P as a good approximation to 
that setup. 



M1-M2 + -7 



1 \ 5 

Ml -M2 + -7j - -7 



(3.11) 



Ml - M2 + -7 > 
Ml - M2 + ^7 < 
Ml - M2 + ^7 = 



3.2 Wavefunctions for the SO(12) point 

We now determine the wavefunctions for the fields localised near an S0(12) enhancement 
point. For that we solve eq. (I2.25P with the above background, closely following the method 

^*The relative signs between the components in expressions p.9p . (|3.10p and p. lip are determined from the 
group theory charges of the states while the overall sign for each curve is determined by studying the form of 
the wavefunctions in section [331 such that given the sign of the fluxes the correct state localises. 
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described in section 12.31 The solution to eqs. (j2.33p can be shown to have the general 
expression 

^ = f {-k2Zi + kiZ2) e-Pil^il'-P2k2l'+P3Si^2+P4J2^i ^ (3_13) 

where km, m = 1,2, and pi, i = 1,2,3,4, are constants depending on the background. In 
this expression / is a general holomorphic function of the particular combination of variables 
that is denoted. In particular, depending on the form of /, the wavefunction can become 
delocalised along one real direction within the matter curve. We return to this point soon 
but first we note that locally we can always take / to be of the form 

fiz) = z' , (3.14) 

where I counts the massless Landau level degeneracy, which is the number of generations and 
so ranges over / = 0, 1, 2. This choice is well- motivated from the existence of local geometric 
selection rules in the evaluation of cubic superpotential couplings, as we discuss in section HI 

Thus, from eq. (j2.36p we have that the ground state wavefunctions can be generically 
expressed as, 

= |P {-k^Zl + klZ2)^ e-PlklP-P2k2p+P35l^2+P452^1 



AT/ 
p 



where is such that 



/ %[^l>l)^ = 5^^^^. (3.16) 

This normalisation factor can be evaluated using standard gaussian formulae, a particularly 
useful such integral takes the form 

I {ni,n2,n3,n4;pi,p2,P3,PA) = / (ww)"'^ (uu)'^'^ {uiv)"-^wu)''^ e 

Js 



«4 p—pi\w\^—p2\u\^+p3UW+p4WU 

(3.17) 



where the rii are positive integers. This integral is simply evaluated to be 

/-A- \ f— l)"l+"2 

I {ni,n2,n3,n^;pi,P2,P3,P4) = T^^RUllidpT' • (3-18) 



\i=l 



PlP2 - PsPi 



In particular. 



vri^g^/l! (IfcipRe(pi) + |fc2pRe(p2) - Re [^1^2(^3+^4)])" ,,, „ 

= — m i+i ll?pll ■ (3-19) 

2 2 (Re(pi)Re(p2)- i|p3+P4P) 2 

An important observation regarding the form of the wavefunction (j3.13p . with / given in 
eq. (I3.14p . is that the combined effect of the Higgs profile and the flux generically induces 
localisation along all directions in S. However, complete localisation does not always occur 
and this can be seen by the fact that the arbitrary holomorphic function when extended over 
the full matter curve can cancel the exponential localisation in one real direction thereby 
delocalising the wavefunction along that direction. The expression (I3.14p is obtained as the 
leading term of a local expansion of the holomorphic prefactor of the wavefunctions. It is still 
valid for use within the triple overlap since the combination of the wavefunctions of different 
matter curves completely localises that integral. However, for the wavefunction normalisation 
this potential partial delocalistion introduces an ambiguity. As we discuss in appendix [Bl for 
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fluxes of order one, as we take in our examples, this ambiguity in the amount of locahsation 
along the curve is not very large. Therefore in the main text we use the normalisation form 
(j3.19p and in appendix [B] we present a more quantitative analysis of the effect on our results 
that a possible wavefunction delocalisation along one real direction can have. 

In what follows we work out the precise expression for km and pi in eq. (j3.15p in terms 
of the background, for each of the curves in the S0(12) enhancement point. Since we have 
already described the procedure with great detail in section [XHl here we are rather schematic. 
Readers which are not interested in the precise form of the parameters may want to skip the 
remaining of this section and jump directly to section HI 



3.2.1 Wavefunctions for the 5m matter curve 



For matter localised on the curve zi = the gauge covariant derivatives appearing in eq. (|2.25p 
are 



Oi = ^ (a + Mi2i) , D2 = ^|S2-Mif2) , 



D:, = -^,,, (3.20) 

Vl 



where Mi = Mi + qyI, and qy denotes the hypercharge of the localised state. Hermitian 
conjugation is obtained in these conventions by complex conjugation in addition to flipping 
the charge of the mode, q ^ —q- 

The matrix B, deflned in eq. (|2.28p . reads 



(3.21) 













'A 







-2Mi 





Vl 










2Mi 







\o 


Ji 

Vl 





0/ 



with dimensionless eigenvalues 



A^^ = 



-pi - Ml , A^^^ = pi - Ml , a;^^ 



2Mi 



(3.22) 



and where pi = \ M^ "I" ( ^ ) • "^^^ corresponding eigenvectors are, 



?0 







/ 





, \ 









-pi - 


Ml 
















V 


K 

111 


/ 



?2 



Pi -Ml 

R 

Vl 



S3 





1 

Vo/ 



(3.23) 



Solving for eqs. ()2.33p then leads to the wavefunctions for the ground state of each of the 
four towers of fields localised in the 5m matter curve, which are given by eq. (I3.15[) with 



Pi = Pi 



P2 



Ml 



1 



5 

P3 



Pa 



(3.24) 



where, following the discussion in the previous subsection, we have assumed Mi to be positive. 
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3.2.2 Wavefunctions for the 10^ matter curve 

The simplest way to deduce the wavefunctions for the 10m curve along Z2 = 0, is to use a 
particular symmetry of the equations of motion (j2.22p for our choice of background. This 
symmetry acts with the following transformations 



Zl O Z2 , 



Ml o M2 , 



Wl -f-)- f2 , 



Di^D2 



(3.25) 



where M2 = M2 + qY'J- Acting with this symmetry on the 5 a/ wavefunctions directly gives 
the wavefunctions for the 10m curve. These are given by eq. (j3.15|) with 



Mo 



P2 



1 



10m 

P3 



Pa 







(3.26) 



and P2 = ■\J M2 + y^j ■ We have assumed M2 to be positive, following again the discussion 
in the previous subsection. 

Dimensionless eigenvalues and their corresponding eigenvectors are given respectively by 



'^0 







-P2 - M2 , 



^2 



P2-M2 



A 



2Mo. 



(3.27) 



and, 



/:10m 
?0 



/A 







( \ 



-p2 - M2 

K 

V2 ' 



I \ 



P2-M2 

K 



1 




(3.28) 



3.2.3 Wavefunctions for the 5h Higgs curve: non-vanishing flux density 

As we have already commented, we should consider three different possibilities for the wave- 
functions localised in the 5h Higgs curve, depending on whether M12 is smaller, bigger or 
equal to zero, with 

M12 = Ml- M2- QYl ■ (3.29) 

In this subsection we compute the wavefunctions for the first two cases, whereas the third 
case is addressed in the next subsection. For concreteness we take M12 < in eqs. (j3.1ip . We 
follow exactly the same procedure as in previous subsections. However, we find convenient 
to express the result in a different coordinate system 



w 



^(^1 + ^2) 



V2 



(zi - Z2) 



V2 



("01 - V's 



In these coordinates the Higgs curve is given by the equation w+u + V-W 
of motion in the (n, t(;)-basis are again given by eq. (j2.25p but 



(3.30) 



0. The equations 



/ 





Du 




D3 \ 




-Du 







-Di 




-D^ 
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(3.31) 
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Gauge covariant derivatives are, 



Dn 



du - M12W 



Rn 



a„ - M 



12U 



D3 = M*i?j 



w u 
— + — 



where we have introduced the following definition 

1 _ 1 ( ^ ± ^ 

v± \/2 V^i ^2 



(3.32) 



(3.33) 



Recall that the orientation is fixed as V2 > vi so that v- > 0. In this basis the matrix 
reads 

/O \ 



Ml 
Rl 







2Mi2 



2Mi2 

with eigenvalues given by the three roots of the cubic equation, 



(3.34) 
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T52 r)2 



12 



■M12 = 



p= 1,2,3 



(3.35) 



and Aq^ =0. The corresponding eigenvectors are, 

/ P 



/1\ 







jJ'H 

(p,w) 



p = l,2,3 



where 



-R 



2M 



12 



k 



5h 

{p,w) 



-R 



2M 



12 



(3.36) 



(3.37) 



With this information at hand we can solve eqs. (j2.33p . For M12 < we find that the ground 
state wavefunctions for each of the four towers of fields localised in the 5// curve are given 
by 

t5, 



P 



N'- 
i\p 



-k^"w + kl"u] e- 



—Pi \w\ \u\ +p^^ wu+p^^ uw 



(3.38) 



with kl« = kl"^^y kl« = kl"^^ and 



5jf 
Pi 



P3 



5h 
P2 



M12, pI' 



{2Muv+V-+Xi"vl)Xi" 
{2Muvl+Xl"v+V-)xl" 



(3.39) 



12 • 



Wavefunctions for heavier modes can be obtained by acting with the gauge covariant 
derivatives ()3.32p on the ground state wavefunctions, as in eq. ()2.37p . The corresponding 
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masses are given in eqs. (j2.38p . Thus, for M12 < there is a massless chiral fermion coming 
from the 5h- Wavefunctions for fields in the 5h representation are also given by the same 
expressions, but their masses are shifted in such a way that there is no zero mode, in agreement 
with the general discussion of subsection 12.31 

Similarly we can work out the wavefunctions for the case on which M12 > 0. In that case, 

XI«XI"XI" = > , (3.40) 

as can be deduced from the determinant of the non-trivial part of eq. (I3.34p . and there 
are one positive and two negative eigenvalues. Thus, wavefunctions are again given by the 
same expressions as in the case M12 < 0, but with the change M12 — s- — M12 and where 
A^^ now denotes the positive eigenvalue. The roles of the 5h and the 5h representations 
is also exchanged in such a way that the massless chiral multiplet now comes from the 5h 
representation. 



3.2.4 Wavefunctions for the 5h Higgs curve: vanishing flux density 

Let us now consider the case on which M12 = for the triplet mode on the 5h curve. In that 
case the eigenvalues of the matrix (i?||/M*)^B in eq. ()3.34p become 



= X^" = 



(3.41) 



and the corresponding eigenvectors 
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(3.42) 

In particular, [Z)3,i)3] = and there is a conserved complex Kaluza-Klein momentum asso- 
ciated to that commutator instead of a quantum harmonic oscillator. We denote by /ckk the 
conserved quantum number. 

Wavefunctions are obtained by following the same procedure as in previous sections. Thus, 
we obtain 

i+n/ f)\\ni 



■p,{n,m,kKK) 



with ground state wavefunction 



(i?||/M,)™+"(DT)"(Z)2) 



VmW.i-Xl")^/'^ (x 



(3.43) 



p,(0,0,fcKK) - Jf^^P 



v1_ + v'^ 
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\v+u + v-wl"^ + 2ilm [kKK{—v-u + v+w)] 



(3.44) 
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where Np is the normahsation constant. The masses of these modes are 



Ml 



Ml 



1. ^kk) 



Mi 



'^2,(n,m,fcj^p^) 



Mi 



'^3,(n,m,fcj^p^) 



(m + n + l)Ai^ + |A:kkP 

(m + n + 2)Af^ + |A:kkP 
(m + n + l)Af^ + |A:kkP 



(3.45) 



Wavefunctions are only locahsed along the directions which are transverse to the Higgs curve 
and arrange in vector-like pairs. Thus, the normalisation constant Np and the correct quanti- 
sation of /ckk are only globally determined. A rough estimate of Np can be given by assuming 
the Higgs curve to be a completely homogeneous space, resulting in 



N„ ~ Rii 



irv-i-V 



-VolK 



1 / TTV^V-R 



1/2 



2R ) e\ 2 

where Volg^ denotes the volume of the Higgs curve and e was defined in eqs. (|2.42|) . 



(3.46) 



4 Wavefunction overlaps 

Our aim is now to compute the coefficient of the cubic couplings in the 4-dimensional effective 
action which results from dimensionally reducing the 8-dimensional theory of section [2Tl The 
relevant operator in the 8-dimensional theory is given by [5] 

VFd / AaAA*. (4.1) 

Js 

Thus, having calculated the wavefunctions for the massless and massive charged fields, cubic 
couplings in the 4-dimensional effective theory are given by integrating the triple overlap of 
the associated wavefunctions over the internal 4-cycle S. Since wavefunctions are localised 
within the local patch so is their overlap, which means that effectively we can perform the 
integral over S as an integral over C^. 

Physical cubic couplings are determined from overlaps of normalised wavefunctions, the 
normalisation condition assuring 4-dimensional canonically normalised kinetic terms. Wave- 
functions for the charged fields are given by 

$12 = x^""'*^" + + X^''*^'' , (4.2) 

where the matrices , t}'^^ and are along the corresponding bi- fundamental generators 
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in su(5) X u(l) X u(l) C so(12). Substituting into (|4.ip we get 



Y 



p,{n,m,l) 



1 

6 











p,{n,m, 


^ 2 p,{n, 



2p,(n,m,l)^l 



X 



10m J 



^ lp,(n,m,l)^ 2 

5ff ^i.5M,i„i,'i-0M,j 



+ 



Off ?/; l/; ' 



(4.3) 

In this expression we have dropped the overall group theory factor Tr | t^^' ,t^^'^' and 
reinstated the generation indices i and j counting the massless Landau level degeneracy and 
the quantum numbers (n, m, I) labeling massive modes on each of the three towers p = 1, 2, 3 
of fields localised in the 5/^ curveEl For the case discussed in subsection 13.2.41 where the 
effective local flux density vanishes, relabeling {n,m,l) — )• (n, m, ^kk) is understood in the 
above expression. 



Note that Y^'^^'^^^ implicitly depends on the hypercharge of the fields involved in the 
coupling. Yukawa couplings correspond to the case of three massless fields and are thus 

given by ^/}o,o,o)' ^^^ere the wavefunctions V'i^(o,o,o)' ^2f,(o,o,o) ^""^ ^i,To,o,o) ^^^^ = 
—1/2. On the other hand, triple couplings between two massless matter fields and one 
heavy coloured triplet are given by l^Jj^^,;), where now ^^1^(0,0,0)' ^ff, (0,0,0) ^nd Xi,|o,o,o) 
have qy = 1/3. Integrating out the three towers of massive (anti-)triplets leads to non- 
renormalisable operators in the infrared, which often put strong experimental bounds on the 
particular F-theory GUT model. Performing this integration is however beyond the scope 
of this work, especially since for the case of proton decay through dimension-five operators 
it would also require knowledge of the up-type cubic couplings localised at Eg enhancement 
points of S. Instead, in this section we analyze the down-type cubic couplings, eq. (14. 3|) . 
between two massless matter fields and the leading contributors in the three towers of massive 
triplets. We split the analysis into two subsections, concerning the cases where the Higgs curve 
has non- vanishing flux density and where it has vanishing flux density and therefore massive 
flelds carry some conserved KK momentum. 



4.1 5// Higgs curve with non- vanishing flux density 

We first consider the case where triplets feel a non-vanishing fiux density and therefore all 
massive modes are localised Landau levels with no conserved KK momentum. The relevant 
wavefunctions were presented in subsection 13.2.31 (and the more general versions in appendix 

In order to evaluate eq. (|4.3p . it is useful to recall some geometric selection rules that 
identify the particular Landau levels coupling to each matter generation. For that, we can 
assign a global U(l) charge to each 4-dimensional field corresponding to the following local 
holomorphic rotation of the internal space 

zi ^ e'^zi , Z2 ^ e*^22 . (4.4) 



^^Cubic couplings between two massless matter fields and the vector multiplets ^o,(n.m,o vanish because of 
the A/" > 1 supersymmetry preserved by the background. They can only be generated after supersymmetry 
breaking. 
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For constant fluxes Mi, M2 and 7, the exponential factor of the wavefunctions is invariant 
under this transformation and so the wavefunctions transform with a phase. The charge of 
the wavefunction is such that in the polynomial prefactor each power of a holomorphic coor- 
dinate zi or Z2, contributes a +1 charge while each power of an anti-holomorphic coordinate 
contributes a —1 charge. 

Massless wavefunctions only involve holomorphic prefactors and so can only have positive 
charges. The charge is essentially labeled by the generation number since the holomorphic 
polynomial prefactors are given by holomorphic coordinates raised to the power of the gen- 
eration index. As studied in [T0l[T2Hl4t [T6] the heaviest generation is usually associated to 
the constant prefactor and so to vanishing charge, whereas second and first generations are 
associated to charges -|-1 and -\-2 respectively!^ 

Massive Landau replicas are obtained by acting with the raising operators on the massless 
fields, as described in subsection l2.31 Ladder operators carry a definite charge under the above 
U(l) symmetry. More precisely, for creation operators we have 

Qu(i) {d\) = +1 , Qf;(i) (^2) = -1 , Qu{i) (^3) = -1 , (4.5) 

with opposite charges for the annihilation operators. Thus, the global U(l) charge of a 
massive field is given bj0 



Qu{i) (^;,{n,m,/) )=i + n-m-l. (4.6) 

The usefulness of assigning such charges becomes apparent when considering the triple 
wavefunction overlap integral. Since the integration measure is invariant under this symmetry, 
the product of the three wavefunctions must also be invariant in order for the integral to be 
non- vanishing. Thus, the only non-vanishing couplings Yp^^ij^^ are such that i+j+n—m—l = 
0. This in particular implies that all Yukawa couplings vanish except for the one of the b 
quark, that is ^/^q^q 0) [inillHllll]. It is however important to emphasise that this U(l) 
symmetry is broken by non-trivial local metric or flux profiles, closed string fluxes [13J or 
non-perturbative effects [42]. We should therefore think of this symmetry as applying at 
leading order in an expansion in the spatial variation of the metric and fluxes. 

The above selection rule also constrains the possible couplings between two massless 
matter flelds and one heavy (anti-)triplet. We have summarised in tabled] the allowed cubic 
couplings between right-handed quarks and squarks and the lightest massive anti-triplets 
'^p,{n,m,i)- The remaining couplings to heavy anti-triplets allowed by the selection rule are 
formed by acting on these lightest states with 'vector-like' combinations of raising operators, 
in the sense of adding no net U(l) charge. Thus, for instance ^r^r couples to the infinite set 
of heavy anti-triplets 

^p,(0,0,0) ; ^p,(l,l,0) > ^p,{l,0,l) ' ^p,{2,2,0) ) ^p,{2,l,l) ; ?p,(2,0,2) ; • • • (4-7) 

and similarly for other pairs of right-handed quarks and squarks. 



^''Another theory of flavour was proposed in [43] for which the generation structure cannot be studied locally 
at an SO (12) point. We discuss the implications for this model in subsections 14.21 and 15.21 

^^Note that, since we only require one generation for the Eh curve, we assume that the global structure is 
such that the generation index vanishes for this curve. 
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UR 


Cr 


tR 




T(OA,0), 7(0,3,1)' 2"(0,2,2); 7(0,1,3), 7(0,0,4) 


7(0,3,0), 7(0,2,1), 7(0,1,2), 7(0,0,3) 


7(0,2,0), 7(0,1,1), 7(0,0,2) 


SR 


7(0,3,0), 7(0,2,1), 7(0,1,2), 7(0,0,3) 


7(0,2,0), 7(0,1,1), 7(0,0,2) 


7(0,1,0), 7(0,0,1) 


bR 


7(0,2,0), 7(0,1,1), 7(0,0,2) 


7(0,1,0), 7(0,0,1) 


7(0,0,0) 



Table 1: Non- vanishing cubic couplings Yp(n mi) between right-handed quarks and squarks 
and the lowest laying massive anti-triplets ?], „j according to the geometric U(l) selection 
rule. We have omitted the subindex p in order to simplify the notation. Similar couplings 
are possible between up-quarks, sleptons and the lowest laying massive anti-triplets. 



In what follows we perform the explicit calculation of the couplings which are shown 
in table [H For large values of R, these are the most relevant couplings, not only because 
they involve the lightest triplets but also because cubic couplings to fields involving a larger 
number of raising operators are suppressed by higher powers of R, as we show below. 

It is worth noting that in the limit where the local U(l) symmetry is exact, and where the 
up-type triplet coupling associated to a point of Eg and the down-type coupling of the S0(12) 
point are completely coincident, so that they obey the same local geometric selection rules, 
we expect that dimension-five proton decay operators do not arise. This is because all the 
massless modes have a positive or vanishing charge, which means that only the dimension- 
five operator involving the heaviest generations can be U(l) neutral (or put another way, the 
coupling between massive up- type and down- type triplets is forbidden by the U(l) selection 
rules). But this operator vanishes by the anti-symmetrisation of the colour indices. The 
relation between the local geometric U(l) symmetries of the S0(12) and Eg points is beyond 
the scope of this work and so we cannot quantify this possibility further. However, we note 
some general things. Firstly it is natural that the local geometric U(l) symmetries of the 
S0(12) and Eg points be strongly correlated as this is required for an appropriate CKM 
matrix. However we also know that the U(l) symmetry must be broken at scales of order 
~ 0.2 in order to induce significant generation mixing in the Yukawa couplings. Therefore 
we do not expect a significant suppression from such a symmetry. Further note that the 
symmetry may be broken even more strongly in the massive sector: although in the case 
studied in this section where the Higgs curve has non-vanishing flux the massive modes are 
Landau-levels and so have a definite U(l) charge, in the next section we study the case where 
the Higgs curve has vanishing flux and then massive modes carry a non-trivial conserved 
KK momentum which completely break the U(l) selection rules. In that case we expect no 
suppression of dimension-five operators due to the local geometric U(l) symmetry. 

Starting from ()4.3p we can pull out an overall factor by expressing the wavefunctions ipi , 
1^2 ™d X ill terms of their respective functions (p, as in eq. (I2.36p . This gives 



y(*.i) ^ 

p,(n,m,l) 



/M,Vio«.vf„"„,) , (4.8) 



where we have defined 

= (5f"^)"(4'")'"(4'")V'" , (4.9) 
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and the overall factor N™^^^ is given by 



n: 



cubic 



1 



P 



R 



P2 +Af2 



R 



+ 4Mi 



12 



2 



Pi + Mi 



Pi + Mi P2 + M2 



(4.10) 



To calculate the couplings of table [T] we rewrite the polynomial prefactor of the wavefunctions 
(^Sa/.j^ (^iOa/j^ V'^rfmi) ™ ^'^^ {Ujw) Coordinate basis (c.f. eq. (13.30p ) and collect powers of 
WW, uu, uw and wu. Making use of the integral (|3.17p we can then express the relevant 
coupling constant as 



Y 



jY-cubic|-_2'jp+g 



r,(0,m,l) 



p q m I 

ka=0 fei, = Ofcc = Ofcd = 
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m/2 



ka I P \ I Q \ I 'm\ / I 



m 



(dl 



1/2 



m—kc 



dl\ (dlA (dl 



Hp + Q- max [ka + kb, k^ + k^] , min [ka + kb, k^ + kd] , 5kQ {5k) , 

-6kQ{-6k) ; pi,p2,p3,p4) 



(4.11) 



where is the Heaviside theta function, 6k ^ ka + kb — kc 
quantities 



Pi = Pi"" + ^ (pi + P2 + Ml + M2 
P3 =pt" + 1 {-Pi + P2 + M1-M2 



p2 = Pa"" + ^ (Pi + P2 + Ml + M2 



kd and we have defined the 
(4.12) 



5h 



1 



p4 = + 2 1 -p^ + P2 + Ml - M2 



as well as 
di 



kfl^)[pl"-Mi2)-Pt"k: 



'^(i,w) [Pi 



12 
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v+ 


pI" 




R 



4Mf2 



(4.13) 



This expression specifies the cubic couplings to the lightest fields localised in the 5h Higgs 
curve. Given some matter generation choice specified by {p, q) the relevant couplings are all 
the combinations of (m, /) such that m + I = p + q, as summarised in table [TJ Although the 
expressions are rather cumbersome they are simple to evaluate with the aid of a computer. 

It is interesting to recall the dependence of eq. (|4.1ip on the local scales iiy and R± or, 
equivalently, on the parameters e and R introduced in eqs. (|2.42p . Such dependence can 
be obtained from the expression of the normalisation constants, raising operators and the 
integral (j3.17p . In particular, in the limit R ^ 1 ^ e where the local approach that we are 
using becomes reliable, eigenvalues scale as in eq. (I2.43P and it is possible to show that 
the above cubic couplings scale as 



Y 



r,{0,m,l) Jl^+l+- 



(4.14) 
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The fact that cubic couphngs to heavy triplets are suppressed by higher powers of R\\R± 
than the Yukawa couphngs has important phenomenological consequences for proton decay, 
as we discuss in detail in subsection 15. ll^i 



4.2 5h Higgs curve with vanishing flux density 

In the case where the effective local density of flux felt by the triplets in the 5h curve 
vanishes, the triplets carry some conserved KK momentum, as discussed in subsection 13.2.41 
Their wavefunction overlap with massless matter fields may therefore differ significantly from 
the case with only Landau levels that we have discussed in the previous subsection. There are 
two key differences: the first is that geometric U(l) selection rules no longer apply since KK 
wavefunctions explicitly break this symmetry. Moreover, the overlap integral now includes a 
new exponential factor depending on the KK momentum. 

The relevant triple overlaps are given by 



Y 



r,(n,,m,A:KK) 



5h 



n/2 



\5h 
■^2 



m/2 



where we define 

r^r,(7i,m,fcKK) 

and the overall factor N^^^^"^ is now given by 



kk) ' 
(4.15) 
(4.16) 
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+{-lY (v: 



cubic 



R 



V2' 



Vl 



. + (pi + Mi 

{V+ - V- ) (p2 + M2 
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V2 



v+) (pi + Mi) 
r = 1,2 

{v^ + v+) (pi + Ml 



(4.17) 



To evaluate the above expression it is convenient to introduce the integral 

/kk ("-1,^2, 71-3,77,4 ; Pi,P2,P3,PA ; 01,02) 

= / ^|;"l'u;"2^"3,j2"4g-Pl|'"'P-P2|«P+P3«'iS+P4W'!«+2jIm(aiio+a2n 

Js 



rU- 



1 



.n2+n4 on 



'-^m ^n-, ^rtr-, 



,714 

ai "^ai 



kll^P2+l°2l^Pl+°l''2P3+''2°lP4 
P1P2-P3P4 



(4.18) 



_PlP2 - PsPi 

In order not to overload this section, here we just present the explicit analytic expression for 



Since the suppression of lighter generations coupling is related to the higher Landau-level number of the 
massive mode it is natural to question how this would be modified by fluxes which do not have a constant 
local profile and therefore break the local geometric U(l) selection rules allowing the lighter generations to 
couple to lower Landau- levels. We expect that the coupling would still be suppressed in a similar, and slightly 
stronger due to the small parameter associated to the flux spatial variation, fashion. Evidence for this can be 
found in section [42] where there is no local U(l) selection rule and we indeed find such a suppression for the 
lighter generations. 
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the particular case on which one of the two quantum numbers of the heavy triplet vanishes 

.{p,g) _ jvcubic(_i)p+g 



Y 



r,(n, O.fcKK) g ,„ E , / E \n/2 



p q n 



ha hj) he 



(-i)Mf)(,!)(:Mter'-fc^'' 



Ikk ^ + q - ka - kb + kc,0, ka + kb + n - kc,0; pi,p2,p3,p4 ; 
where have introduced the quantities 



kKKV+ -kKKV- 



2RJvl + vl 2RJvl + 

and pfc are given as in eqs. (|4.12|) with the replacements 

p'"^ "i"^ (420) 



Explicit analytic expressions for the couplings to other heavy triplets can be worked out in a 
similar fashion starting from eq. ()4.15p . 

An important qualitative feature of the coupling constants ()4.18p is that there is an ad- 
ditional exponential suppression with respect to the case with no KK momentum discussed 
in the previous subsection. This suppression, however, cannot be significant for the lightest 
KK modes if we are to keep within the local approximation, as studied in section 12.41 It 
is simple to check this explicitly by varying the parameters. Nevertheless, a more intuitive 
understanding is as follows. The exponential suppression comes from the oscillations of the 
5h wavefunction with KK momentum. If the gaussian decay envelope of the matter wave- 
functions is much larger than the oscillation frequency of the KK state then the oscillations 
cancel, leading to an exponentially suppressed overlap integral. However, within the local 
approximation the gaussian width must be less than the length of the 5h curve and so the 
lightest KK state should have at most one oscillation within the gaussian envelope. Thus, 
for the range of validity of the local approach, this extra exponential suppression on top of 
the polynomial one of eq. (j4.14p is rather mild. 



5 Phenomenological implications 

We now address a more quantitative analysis of the cubic couplings that we have presented in 
the previous section and discuss their possible phenomenological implications. The primary 



^^The expression for Y^^^^^^ kKn) ^'^^^'^^^ from this one by making the replacements d^^ — >■ d^^ and 
using the function 

r I , , 111,1 kKKV+ fcKK«- ) 

iKK I p + q - ka — Kb, kc, ka + kt,n - kc', pi,p2,p3,p4 ; 
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application of our results is to proton decay induced by dimension-five operators. More pre- 
cisely we have computed the coupling of the right-handed quarks and squarks (or equivalently 
of the up-quarks and sleptons) to the massive down-type coloured triplets which mediate pro- 
ton decay. Although this coupling does not contain the full information needed to calculate 
the coefficient of the dimension-five operator, it plays a key role in such a calculation. Sup- 
pressing this coupling is a sufficient but not necessary condition to avoid the present strong 
experimental bounds on such dimension-five operators. One of the most important attributes 
of the coupling is that it is a superpotential term, which means that due to its holomorphic 
nature it cannot involve string oscillators. The latter therefore can only enter in the physical 
coupling through the normalisation of the fields. 

Besides the most direct application to proton stability, the coupling of matter fields to 
heavy modes is also important in the context of other non-renormalisable operators of phe- 
nomenological interest. For instance, in F-theory implementations of the Froggatt-Nielsen 
mechanism [33] Yukawa couplings arise from higher dimension superpotential couplings. Sim- 
ilarly, FCNC operators can often be induced by integrating out heavy fields and may play an 
important role in studies of supersymmetry breaking. We further discuss the implications of 
our results for these non-renormalisable operators in subsections 15.21 and 15.31 respectively. 

5.1 Coupling to heavy modes: proton decay 

As discussed in section [3l the nature of heavy triplets differs according to whether or not they 
feel a non-vanishing flux density on the curve. We consider these two possibilities in subsec- 
tions [5TLT] and [5TL21 In both cases the strategy is similar: we evaluate the leading couplings 
of interest and study their behaviour as we vary the input parameters. Although there are 
many input parameters that go into the calculation, we find that the qualitative behaviour 
depends primarily on the ratio of the fluxes to the parameter R defined in eqs. (I2.42p . Hence, 
a simplified strategy that we adopt here is to fix the flux parameters to some 0(1) values 
and to study the behaviour of the coupling as we vary R for some fixed e < 1. The scanning 
range for R is set by the limits for which the effective theory is under control, as studied in 
section 12.41 keeping in mind their approximate nature. For R < 1 the local approach that 
we have adopted in our computations breaks down and higher derivative corrections to the 
8-dimensional effective action (in the type IIB language, corrections from string oscillators 
to the normalisation of the fields) become important. On the other hand, too large values of 
R are in tension with the observed value of acuT , requiring values of e close to the unity, for 
which higher derivative corrections become again relevant. Thus, given some region of values 
for R, the value of e is fixed within a range by the validity of the 8-dimensional effective 
theory and the relation between the local scales and ogut in that particular model. 

5.1.1 Coupling for Higgs curve with non-vanishing flux 

The first case that we model is where the flux density on the 5h Higgs curve is such that there 
is a massless anti-triplet. Without any other extra ingredient, this scenario therefore does not 
fully account for doublet-triplet splitting. The massless anti-triplet would have to be lifted 
by some deformation of the geometry, which in the effective field theory would correspond to 
a supersymmetric ;U-term or to a vev for a GUT singlet. If the resulting mass is well below 
the local KK scale we expect this to only induce a small correction to our calculations. The 
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lifted massless anti-triplet, however, leads to a new extra contribution to proton decay which 
can in particular dominate over the one of the massive modes that we are computing. In this 
sense the contributions that we consider in this subsection should be regarded as only part 
of the total one in a complete mo del 1^ 

The relevant coupling was calculated in section 14.11 As discussed in the introduction 
our primary interest is in the ratios of the heaviest generation Yukawa coupling to the cubic 
couplings to heavy anti-triplets. This measures the departure from the naive 4-dimensional 
equality of Yukawa couplings and cubic couplings to heavy triplets. In order to present the 
results most concisely we define an average coupling to the dominant massive modes for each 
pair of generations 

i^m) = l+p + q Yl K^!{n!m,i) ' (5-1) 

n,m,l 

where the sum runs over the dominant massive anti-triplets in the r-th tower of localised 
heavy fields in the 5h Higgs curve. See table [T] for an explicit list of the terms. 

In figure [2] we plot the average coupling to the first tower of massive anti-triplets, the 
one generated by acting with raising operators on the massless modes, for values of the flux 
of order one. Since doublet-triplet splitting is not fully accounted for in this case we turn 
off the hypercharge flux for simplicity 7 = 0. The plot shows the coupling to the leading 
massive modes for different generations. The dashed line shows the triple coupling of massless 
modes, i.e. the Yukawa coupling of the heaviest generation. The couplings to the other two 
towers, which do not have a zero mode, are similar in nature and are shown in figure Isf^ For 
completeness we also plot in figure [5] the masses of the different coloured anti-triplets that 
participate in the average couplings (15. ip . 

The behaviour of the cubic couplings to massive modes is relatively universal and robust 
with respect to small variations of the flux and R. We find that for small values of R, such 
that it is of order the flux parameters, the coupling to massive anti-triplets for the lighter 
generations is only slightly suppressed compared to the Yukawa coupling. This means that 
for this region of the parameter space the constraints on proton decay are stronger than the 
naive 4-dimensional field theory result using ()1.3p . As we increase R relative to the flux the 
lighter generation couplings decrease more steeply than the heavier generations, such that for 
large enough values of R the couplings of the lighter generations are significantly suppressed 
compared to the minimal 4-dimensional field theory guess. However, note that the bottom 
Yukawa coupling also decreases for large R leading to very small tan /3 which eventually 
becomes incompatible with gauge coupling unification, which favours tan /3 > 1. Therefore 
in such a setup we expect that from the Yukawa coupling not being too suppressed that 
smaller values of R are favoured. 

^'^Nevertheless, the wavefunctions and overlaps for this setup are identical to an analogous setup where this 
problem is avoided. This is where the hypercharge flux is such that there is a massless triplet rather than 
anti-triplet. Still some deformation of the geometry is required in order to lift the massless mode, but now 
since it is a triplet rather than anti-triplet it does not couple to the massless matter fields through down-type 
Yukawa couplings and does not enhance the rate of proton decay. Therefore we regard this calculation as still 
useful for realistic examples. 

^^The dip seen in the couplings to the third tower shown in figure [3] goes all the way to zero, although due 
to the resolution of the plot it is not completely captured. The origin of this dip is on the overall factor H4.10|l 
which develops a zero as a function of R for that particular tower. A similar feature occurs in figure [B] 
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Figure 2: Average couplings of different generations to massive anti-triplets in the first tower 
of localised fields, for flux values vi = 5/6, V2 = 5/4, Mi = 1.6, M2 = 2, 7 = and e = 1/10. 
For reference, we also show the Yukawa coupling for the b quark (dashed line). 




Figure 3: Average couplings of different generations to massive anti-triplets in the second 
(left) and third (right) towers of localised modes, for flux values vi = 5/6, V2 = 5/4, Mi = 1.6, 
M2 = 2, 7 = and e = 1/10. 



5.1.2 Coupling for Higgs curve with vanishing flux 

We study now the case where triplets in the 5h Higgs curve feel a vanishing total effective 
flux and therefore, given the appropriate topology, there are no massless triplets and doublet- 
triplet splitting is accounted for. Massive vector-like pairs of triplets carry a conserved KK 
charge and their wavefunctions are delocalised along the Higgs matter curve. This means 
that we cannot calculate the wavefunction normalisation explicitly and have to estimate it, 
as we did in eq. (I3.46p . However, this is not a particularly serious shortcoming given the 
accuracy at which we are working, and we refer to appendix B for a more detailed study of 
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Figure 4: Masses of the anti-triplets participating in the average couphngs (15.11) . as sum- 
marised in table [H for the first tower of modes with flux values vi = 5/6, V2 = 5/4, Mi = 1.6, 
Ma = 2, 7 = and e = 1/10. 



this issue. 

We plot in figure [5] the coupling of the different matter generations to the lightest massive 
vector-like pair of triplets in the first tower of localised fields, for values of the flux of order 
one. The flux is chosen such that there is no total effective flux for the triplets, so that M12 
as defined in ()3.29p vanishes. The massless chiral spectrum is determined by the expressions 
(I3.9p - (]3.1ip and is such that that there is a massless chiral Higgs doublet as well as the 
appropriate chiral matter fields. The massive vector- like pair of triplets that we consider is 
the lowest KK mode with vanishing Landau-level quantum numbers. As discussed in section 
[3l the presence or absence of a vector-like pair of massless triplets is only determined by the 
global geometry. From the local approach we assume that such a massless pair is absent and 
therefore all triplets carry non-zero KK momentum. The coupling to the lightest vector-like 
pair of triplets in the other two towers of localised fields is plotted in figure [6l 

The general behaviour of the couplings is similar to the setup discussed in the previous 
subsection. The more detailed differences being that the Yukawa coupling is slightly less 
sensitive to the values of R, allowing for a larger viable range. Also the suppression of the 
lighter generations is stronger than that in figure El Both of these changes improve the 
prospects for suppressing proton decay operators. 

5.1.3 General overview 

Having studied both the possible cases for the Higgs triplets in the last two subsections we 
can attempt to draw some general conclusions from our analysis of the coupling of the heavy 
triplets to the different generations. Of course all conclusions are qualified with the fact that 
our analysis can only be taken to hold up to order 1 factors or so. Also the possible parameter 
space is restricted by our use of a local 8-dimensional field theory as discussed in section 12.41 
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Figure 5: Couplings to the lightest massive vector-like pair of triplets in the first tower of 
fields localised in the Higgs curve, for flux values vi = 5/6, V2 = 5/4, Mi = 1.6, M2 = 2, 
7 = 1.2, e = 1/10, /jkk = 1 and trivial Landau-level quantum numbers. For reference, we also 
show the Yukawa coupling for the b quark (dashed line). Note that, since hypercharge flux 
is non-vanishing the coupling is not GUT group universal. The coupling to the triplets that 
we plot is (3, 1)1/3 (1' 2)_i/2 (3, 2)i/g while the Yukawa coupling plotted for comparison 
is (l,2)_i/2 0(3, 1)1/3 0(3, 2)1/6. 




R 



Figure 6: Couplings to the lightest massive vector-like pair of triplets in the second (left) and 
third (right) towers of fields localised in the Higgs curve, for flux values fi = 5/6, ^2 = 5/4, 
Ml = 1.6, M2 = 2, 7 = 1.2, e = 1/10, /ckk = 1 and trivial Landau-level quantum numbers. 



and it may be that leaving this framework could change the results in a quantitative and 
perhaps even qualitative way. With this in mind however it is also important to note that 
the behaviour we observe, up to order 1 factors, is quite robust against variations of the flux 
parameters and geometric scales within their respective allowed regimes that allow for the 
perturbative field theory analysis we have performed. 
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Figure 7: Proton decay superpotential couplings involving only one lightest generation of 
quarks. The diagrams involving these operators are expected to dominate in the bulk of 
parameter space where the coupling to the heavier generations are only mildly suppressed. 

Perhaps the most general and important property of the couplings we have calculated 
is that there is a suppression of the coupling of lighter generations to the heavy triplets 
in analogy with the minimal 4-dimensional field theory GUTs behaviour, and which for 
large enough values of R can lead to substantial suppressions of the couplings. Another 
general pattern we find is that the bottom quark Yukawa coupling is small, implying a small 
tan /3 regime is perhaps more natural. Both of these properties are favourable in terms 
of suppressing dimension-five proton decay. The suppression for the lighter generations is 
particularly important, allowing to generate the small numbers needed to qualitatively match 
dimension-five proton stability constraints (assuming that a similar suppression for up-type 
couplings holds). 

A quantitatively precise examination of these results would require the precise knowl- 
edge of the model dependent relation between the local scales R± and and the global 
compactification scales parameterised by Mpianck and acuT- For nearly homogeneous com- 
pactifications, where small values of R are most natural, such analysis reveals that within 
the bulk of the allowed local parameter space it is quite difficult to suppress the coupling 
to triplets sufficiently to weaken the proton stability constraints of minimal 4-dimensional 
GUTs that result from eqs. ()1.3p . The strongest constraints on the parameter space come 
from the relation to oqut and from the requirement of a sufficiently large bottom Yukawa 
to be compatible with tan /3 > 1. Both of these are approximate and if we allow ourselves 
to go to the edge of parameter space, say with R ~ 25, then further suppression is possible. 
This suppression of the couplings relative to the expected 4-dimensional values is most sig- 
nificant for the lighter generations which means that the strongest constraints would come 
from superpotential couplings involving as many heavier generations as possible, such as the 
presented in figure [TI^ We cannot rule out that allowing for order 1 factors and going to 
the edge of allowed parameter space, as well as assuming a similar suppression to the up- 
type Higgs triplet couplings, all the proton decay diagrams could be suppressed sufficiently to 
avoid experimental constraints in this class of models. However, while keeping this possibility 
in mind, within the framework we have used the most natural conclusion is that in F-theory 
models that are based on nearly homogeneous manifolds the coupling of the heavy triplets 
is of the same order as, or even more enhanced than, in minimal 4-dimensional field the- 
ory GUTs and further suppression of proton decay operators is required through additional 

^^This information is particularly non-trivial in the presence of additional selection rules which may, for 
example, forbid only the operators involving heavier generations but not lighter ones. 



35 



mechanisms. 



5.2 Relation to Yukawa couplings 

We have computed the couphng of massless chiral modes to locahsed massive modes using 
wavefunction overlaps. So far our primary application of this calculation has been to pro- 
ton decay through dimension-five operators. However the couplings we have calculated are 
important also for other physical quantities. In this section we present a brief discussion of 
their application to Yukawa couplings. 

The nature of Yukawa couplings in F-theory GUTs depends on how the multiplicity of 
the generations arises. So far we have considered the case where the three matter generations 
arise from the multiplicity of Landau-level zero modes of a single curve. Yukawa couplings 
are obtained from zero-mode wavefunction overlaps and, as we have seen, for this setup they 
have a rank one structure due to the local U(l) symmetries [THj. Yukawa couplings for the 
second and third generations should then arise from deformations of this structure. It was 
shown in |12H14j that such deformations cannot arise from worldvolume gauge fluxes, which 
are the only ingredients that we have included so far in our analysis. Instead, to generate a 
higher rank structure locally what is required is a non-commutative deformation of the local 
gauge theory. This can arise from two sources: either through Imaginary-Anti-Self-Dual 
(lASD) closed string fluxes of type (2,1) [13], or through non-perturbative effects such as 
instantons or gaugino condensation on some distant brane [52] . In [65j it was argued that in 
fact the two possibilities are equivalent when it comes to their effect on the local theory since 
non-perturbative effects source lASD fluxes. 

Calculations of Yukawa couplings within this framework were performed in |131ll6j. In 
particular it was shown in |16j that the deformed Yukawa couplings can be written as 

i^totai = Y}^i + 6 + y w ) + ... . (5.2) 

Here the subscripts on the Yukawa coupling contributions denote the type of operator in 
the gauge theory, with 'tree' denoting the usual Yukawa coupling arising from F A ^ as 
studied in this paper, and 'n.p.' denoting the higher dimension operator induced by the 
non-perturbative effect (which takes the schematic form AF). The parameter e is related 
to the non- commutative deformation or, equivalently, to the non-perturbative scale. The 
superscripts on the operators denote to which order in the deformation of the wavefunctions 
the operators are evaluated at. Thus, a superscript denotes the undeformed wavefunctions, 
while a superscript 1 denotes that one of the three wavefunctions in the triple overlaps is 
the first order deformed one. It was further shown that the deformed wavefunction can be 
written as an explicit linear combination of the massive wavefunctions of the undeformed 
theory. These wavefunctions are exactly the wavefunctions we have been calculating and 
the resulting coupling l^^ce directly given by the appropriate linear combination of the 
couplings of two massless and one massive modes which we have calculated. Therefore the 
calculations of massless and massive wavefunction overlaps presented play an important role 
in calculating also the Yukawa couplings. There are two direct applications of our results. The 
first is to use within a Yukawa coupling calculation analogous to [16] but for the more realistic 
S0(12) rather than U(3) setting (though in practice wavefunctions of both setups turn out 
to be rather similar). Secondly we see that there is a connection between Yukawa couplings 
and proton decay. In particular, within the region of parameter space where proton decay 
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operators are suppressed also this contribution to the Yukawa couphng would be suppressed. 
Yukawa couplings would then primarily have to arise from the last term in ()5.2p F^ 

This model of flavour that we have been using is elegant and has attractive features. Be- 
sides its simplicity, it provides a very natural explanation for the observed hierarchy between 
the Yukawa couplings of the third and the two lighter generations. It suffers, however, from a 
number of difficulties in attempting a realistic implementation. In particular, no global com- 
pact stable backgrounds are known, either in IIB or F-theory, which support the required 
non-commutative deformation. Moreover, in order to obtain a realistic flavour structure 
the Yukawa deformation should arise at order ~ 0.2, which seems too large to have a non- 
perturbative origin while still maintaining a perturbative expansion. Our results showing that 
within the local perturbative regime R > 1 wavefunction overlaps are generally suppressed 
are likely to make this problem even more severely 

An alternative model of flavour was proposed in 03] where each generation arises on a 
different matter curve and the flavour structure comes from embedding the SU(5) GUT into 
Eg through the Froggatt-Nielsen mechanism. The calculations that we have performed in the 
previous sections are also relevant for these models, as Yukawa couplings in this case arise 
from higher dimension operators which, in turn, come from integrating out heavy modes. For 
example, starting from a bottom-type coupling to heavy modes and integrating out those, 
generates the s quark Yukawa coupling. The relevant operators are 

W D 5H,5fH0t + X5,5^K + MKK^fHf^ , (5.3) 

where we have denoted the heavy modes with the superscript 'KK'. Integrating out the heavy 
modes gives 

^5„,5,10, , (5.4) 

where X is a GUT singlet that develops a vev. The first coupling of ()5.3p is of the type we 
have been studying, with the small difference that the massive mode is taken for the matter 
curve rather than the Higgs curve. Since the flavour structure in the models of |43) comes 
from the additional U(l) symmetries, it would require a large coupling to the massive modes 
and so in this case the, more natural, small R region of parameter space is preferredEl The 
second term in eq. (j5.3p is the type of coupling that we study in the next section. Finally note 
that in the context of proton decay the singlet vevs already account for the 4-dimensional 
coupling suppression as in (jl.3p . and so only a further mild suppression, coming for the 
difference between the doublet and triplet coupling as we have been studying, is required 
to avoid proton decay constraints. Given our results, which showed that couplings to heavy 
modes can be easily and generically much smaller than 1, such a possibility seems quite 



■^■^ Actually, as noticed in [16], the coefficient of the last term in (|5.2p vanishes for an SO(12) enhancement 
point and the non-trivial contribution to the Yukawa couplings is generated at higher order in the deformation. 
In that sense the suppression of the first term may be quite important. 

There is also a potential difficulty discussed in |15] which is that the initial rank one structure applies 
only locally to a single enhancement point of say SO(12). However, in concrete compactifications there are 
many such enhancement points, which are spatially separated, generically leading to a higher rank structure 
and destroying the hierarchical nature. The presence of a number of enhancement points however is not a 
no-go theorem as not all say SO (12) points must correspond to the down- type Yukawa but can amount to 
intersections of different curves. 

^^Note that in the limit of constant fluxes, due to the local geometric U(l) selection rules, the leading such 
coupling would involve one of the massive towers which does not have a massless mode. 
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natural. However an explicit computation of the dependence on the generation structure 
would require going beyond the local S0(12) framework we have been using. 



5.3 Soft masses and FCNC 

As another application of the type of computations we have performed we study a set of 
operators which may be responsible for soft masses after supersymmetry breaking. This is 
a particularly interesting set of operators to investigate because a well-known criticism of 
gravity mediated supersymmetry breaking is that FCNCs are expected to be generated after 
integrating out heavy modes. The coupling of the massless matter sector to heavy modes is 
therefore crucial in calculating this effect o 

More precisely, let X = . . . + 9^Fx be a supersymmetry breaking superfield and Qi a 
chiral matter superfield (baryonic or leptonic) belonging to the i-th generation. If the Kahler 
potential contains dimension-six operators of the type 



K D K.jxx J d^d X^XQ]Qi , (5.5) 

then soft masses are generated 

rh% = K^,xx \Fx\^ ■ (5.6) 

Denoting by Vq^ij — Sij + eij the unitary matrix which diagonalises the fermion mass matrix 
(with Eij << 1 and ejj = for i = j), the scalar soft masses in the basis where fermion masses 
and gaugino-fermion-scalar matrices are diagonal are 

mjj ~ m| + eikmlj - mleij . (5.7) 

Experimental results on flavour violation in the Standard Model set severe bounds on the 
off-diagonal terms in the family space m?- |67] . 

There are various possible sources for the dimension-six operators ()5.5p within the F- 
theory GUTs context. In figure[8]we present the leading such contributions. The first type of 



Figure 8: Feynman diagrams leading to dimension-six operators of the type (j5.5p by exchange 
of massive U(l) gauge bosons (left) or massive scalars (right). 



contribution, introduced in [49], is obtained by integrating out a massive U(l) gauge boson. 
This contribution has been argued to be flavour universal. However, besides this contribution 







^^For non-perturbatively generated FCNC in string theory see [66| 



38 



there can also be non-universal contributions arising from the exchange of heavy Landau or 
KK states which we study in what follows. 

The relevant cubic coupling constants \jn involved in the diagram of figure [5] (right) are 
given by the overlap integral 

f ^-x^-j^n = {n\XQ,) . (5.8) 
Js 

In this expression we have used a convenient Dirac notation, where actually {x\n) = de- 
notes the complex conjugate of the wavefunction Hence, the coefficient of the dimension- 
six operator eq. (jS.Sp which results from these contributing diagrams is 

„ _ ^mAjn _ sr^ {XQi\n){n\XQj) 

with Mn the mass of the state Such a sum over mediating states has a lower bound 
obtained by keeping the lightest state with mass Mi, 



> ^ • (5.10) 



Similarly, an upper bound can be obtained by using the inequality 

_^ {XQ,\n){n\XQ,) {XQ,\XQj) _ 1 /" ,2^ ^ 

which results from making use of M„ > Mi and the completeness of Landau levels. It is 
important to recall here that eqs. (I5.10p and ()5.1ip put bounds on the particular contribution 
of the Landau states $n to the operator ()5.5p . and not on the operator itself, which in general 
may receive other contributions (for instance from string oscillators) not accounted for by 
these expressions. 

In order to perform a more precise estimation of the quantities appearing in eqs. (j5.10p 
and (15. lip we have also to specify the origin of the SUSY breaking field X. If it is a modulus- 
like field, singlet under the Abelian symmetries we are considering, local considerations do 
not determine its wavefunction and we cannot say much by using the methods of this paper. 
However, if X is a charged field then it actually experiences the Higgs vev and/or the Abelian 
fluxes and it has some degree of localisation. This is the case we consider in what follows. 

Since the cubic coupling we are interested in involves two fields charged under the GUT 
group and one GUT singlet, the S0(12) model that we have been considering in previous 
sections is not accurate for this case. A more suitable toy model based on the gauge group 
U(3) and containing all the necessary ingredients was studied in |16] . that we consider in 
what follows. 

In this model massive state J^'^l are actually labeled by three quantum numbers |^n) — ^ 
l^(m,n,/))- Their masses are given by 

^(L,n,0 = ^(mA+ + n|M|+/|A-|), (5.12) 



^^We identify our fields X, Qi and "l>n with the fields in 16 according to 
M < 0. 
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where = ^ (^M ± ^ M"^ + 4.E?m^^ , M and m^f) are respectively flux and Higgs parameters 

and we have made explicit the dependence on the local scales R\\ and R±_ introduced in section 
[2j We consider the minimal case in which the multiplicity of X and is one, whereas the 
multiplicity of Qi is taken to be three, i = 0, 1, 2 with z = labeling the heaviest generation. 

We denote by A*^^^;^ the cubic coupling between X, ^(^rn,n,i) Qi. As it was the 
case in the S0(12) setup, the local geometric U(l) charge conservation (j4.4p only allows for 
couplings which satisfy i — m — n + l = 0. Thus, the lightest fields coupling to X and to 
different generations of squarks and sleptons are 

1st generation: $(2,0,0) , ^{0,2,0) , ^(1,1,0) , • • • 
2nd generation: $(i,o,o) > ^{0,1,0) > • • • 
3rd generation: *l'(o,o,o) j ■ ■ ■ 

In the toy model that we are considering <^*(o,o,o) is massless. However, we expect that in a 
more realistic setup there would not be such massless exotic state. This mass lifting is not 
under control locally and therefore we have nothing to say about it here. 

For simplicity we choose the fiuxes of the model to be small enough 2M2 < R'^m^, such 
that the states with m = I = are lighter than the other states. In that case cubic couplings 
are of the form 



A^ 

^(0,n,0) - g 



+(VfV?" -^-f"V'f)# (5.13) 



and, in particular, the leading couplings to massive states for the first two generations are 
"^^0 2 0) ^lo 1 0) • Making use of the explicit expression of the wavefunctions in [16] we find 



^1 _ \M\'/'X'{X--V2Rml) I .-|,,|.+.-|.,|._^|.,_.,|. 



2 

|2 I x-U„|2_"Vu,_,„|2 



S 



_ \M\X (A -V2Rml) f A-|.,|2+A-|..P-!||.,-.,|2 
^(°'^'°) - V2N^Ni,NxR^ml is"'"'" 

where N^, Nq and Nx are wavefunction normalisation constants. Using the integral ()3.18p . 
together with eqs. (|5.1U|) and (|5.12|) . the corresponding terms in the Kahler potential are then 
readily found to be 

K_ _ ^ 1^(0.1.0)1' 

'''''' ^(0,1,0) 2MHmNl^NxyR^mlX^X-V2Rml)^ ' 

iA?o,2,o)P _ -'^'m 



^22^^ ' MJ 2,0) 4M2(A-)4(iV^iV2Ar^)2(A- _ V2iim|)4 

Hence, due to conservation of the geometric U(l) charge no off-diagonal terms are generated 
in the Kahler potential. We know however from eq. (j5.7p that this is not enough to suppress 
FCNCs, since unequal terms on the diagonal do generate FCNCs after diagonalising the 
fermion mass matrices. 
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Indeed, from eqs. (|5.15p we obtain 



- 2 ~ 2 
"^22 - "^11 



mlM^{X'f{\--V2RmlY 
where we have used the normahsation of the wavefunctions 



1/2 



\M\j R'^mlM 



\Fx\' , 



(5.16) 



(5.17) 



For a suppression of the FCNCs the two mass terms mf^ "^22 would have to be equal to 
a high accuracy. 

A similar reasoning can be followed for the upper bound (j5.1ip on the Kahler term. 
Starting from the cubic couplings (|5.13|) and integrating out over the massive states, one has 



K - < ?<'™a,x_ 
J^ijXX ^ ^ijXX 



which, making use the explicit form of the wavefunctions, becomes 



^^(0,j,0)^(0,i,0) Js 



(5.18) 



7>'max_ 
^ijXX 



3(A- 



V2X- 



2R'^mi 



+ 1 



zlz{ e 



Rml 



(5.19) 



This integral can be computed as before. We do not write here its explicit analytic expression. 
Notice however that, similarly to the lower bound that we have just discussed, the result- 
ing Kahler terms are again diagonal in the generation space due to geometric U(l) charge 
conservation, but, importantly, not degenerate. 

Since this is a toy model we do not perform an analysis of eqs. (l5.15p and ()5.19p on their 
parameter space. Our main emphasis is to show how contributions of the form (j5.5p can 
be calculated and that indeed for generic parameters they are not flavour-universal and lead 
to large FCNCs. In particular, there is a priori no natural reason why the flavour-universal 
contributions from U(l) gauge boson exchange should dominate the contributions studied 
here, since the masses of the corresponding gauge bosons and of the scalar fields are expected 
to be at the same scale. However, having shown how such contributions can be explicitly 
calculated, it would be very interesting to explore the full parameter space, particularly 
within a more realistic model, and see if there are regions where the FNCN operators can be 
sufficiently suppressed. 

Finally, let us raise a further point regarding the U(l) boson exchange contribution itself. 
The flavour universality in this case arises because the gauge boson is taken to couple in the 
Kahler potential in such way that diagonalising the kinetic terms for the matter fields also 
diagonalises the coupling to the gauge boson. However, in general this is only possible if 
the wavefunction profile of the gauge boson is constant. We therefore expect flavour non- 
universal couplings to higher KK states of the U(l) gauge boson, as the wavefunctions of 
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these states are not constant. Once more studying the nature of these FCNCs reduces to 
studying massive mode wavefunction integrals, but this time for the U(l) gauge bosons. Such 
studies are however beyond the scope of this work. 

6 Summary 

In this work we have studied the couphng of matter fields to heavy modes through cubic 
superpotential interactions in the context of SU(5) F-theory GUTs. We have considered, in 
particular, fields that are localised within a patch around a point of S0(12) enhancement and 
computed the local form of their wavefunctions, for both massless and massive modes. Down- 
type Yukawa couplings are obtained by integrating overlaps of three massless wavefunctions, 
whereas couplings of two matter fields to the heavy coloured triplets are given by integrating 
overlaps of one massive and two massless wavefunctions. 

The coupling of MSSM fields to heavy triplets is an important ingredient in the study 
of proton stability through dimension-five effective operators that result from integrating 
out the heavy triplets. In this context our calculations are a prerequisite to constraining F- 
theory GUTs through dimension- five proton stability. We find that, analogously to minimal 
4-dimensional GUTs, the coupling of the heaviest generation to the massive triplets is of the 
same order as the Yukawa coupling, and that the coupling of the lighter generations to the 
triplets are suppressed with respect to the heavier ones. The quantitative analysis depends 
on few local parameters associated to the flux and geometry, which are in turn constrained 
by the measured values of acuT and Mpianck and also by keeping the bottom quark Yukawa 
not too small so that we can remain within the tan /3 > 1 regime favoured by unification. 
The detailed relation between the local scales and the global ones is however very complicate 
and model dependent, and given a particular compactification it is in general out of the scope 
of the present techniques to obtain such relation in a precise way. Once the relation between 
local and global scales is understood in a given model, the constraints on the local parameter 
space coming from the observed values of Mpianck and ogut would be made precise within 
that model. 

For the particular case of nearly homogeneous compactifications we find that the relation 
between local and global scales is such that for the most natural values of the local parameters 
the couplings to the heavy triplets is typically the same or, in the case of the heaviest gen- 
erations, stronger than in minimal 4-dimensional GUTs. Thus, constraints on the coefficient 
of dimension-five proton decay operators involving heavy generations in this class of models 
are stronger than in minimal 4-dimensional GUTs. We can also identify regions on the edge 
of local parameter space where the couplings to the triplets are more suppressed than in 4- 
dimensional GUTs for the lighter generations, thus weakening constraints from proton decay 
on those operators. However, within these region tan /3 is expected to be very small. 

The computation that we have performed is only a part of the complete computation for 
dimension-five proton stability since cubic couplings to up-type triplets localised at points of 
Eg enhancement also participate in the above operators. We expect a similar computation of 
those couplings to be technically more involved than the one that we have carried out here, 
mainly because of the presence of local monodromies. It would certainly be very interesting 
to perform such a computation. In particular, note that the up-type ratios of couplings anal- 
ogous to eqs. ()1.3p are actually significantly smaller than the down- type ones, and therefore 
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even stronger suppression in the F-theory couplings near an Eg point would be required to 
match or to overcome the suppression in minimal 4-dimensional field theory GUTs. 

The coupling of massive modes to massless ones has a number of phenomenologically 
interesting applications, some of which we have already discussed. Within flavour physics 
the computation features in that of Yukawa couplings as studied in [16J. It is also a key 
computation in realising the Froggatt-Nielsen model of flavour proposed in [43j, where the 
higher dimension operators come from coupling of two massless modes to one heavy mode. 
Interestingly, since in the above model of flavour the suppression of the coupling to the heavy 
Higgs triplets is due to the U(l) symmetries at least as strong as in minimal 4-dimensional 
GUTs, only a slight suppression of the coupling to the triplets relative to the doublets is 
required to evade constraints from proton decay. Given that the coupling to heavy triplets 
can be naturally suppressed it would be interesting to study if the required small suppression 
could be realised within this framework. 

Another application of this type of couplings that we have discussed is in understanding 
FCNCs induced by supersymmetry breaking. Within a toy model we were able to give 
expressions for the magnitude of FCNCs induced by integrating out heavy modes. We flnd 
that FCNCs are induced by tree-level diagrams and are suppressed by the mass scale of the 
heavy modes. In particular we find no natural suppression relative to the flavour-universal 
contribution coming from heavy U(l) gauge bosons studied in |49) . 

There are further phenomenological applications that one might think of, such as the gen- 
eration of neutrino masses through the seesaw mechanism. Since the right-handed neutrinos 
are GUT singlets the neutrino Dirac masses would be of the same form as the coupling studied 
in section 15. More generally, we believe that the local approach opens up a calculational 
framework for studying the interactions with massive modes within a phenomenologically re- 
alistic and very general setting. This is important because it is precisely in understanding the 
interaction between infrared and ultraviolet physics that string phenomenology as a subject 
has its strongest claim as a necessary tool in physics. 
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For the type I seesaw mechanism the Majorana nature of right-handed neutrinos would however require 
different techniques to the ones that we have been using in this work. 
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A Wavefunctions for oblique hyper charge flux 



In this appendix we consider a slightly more general set of U(l) fluxes, which in addition to 
eq. ()3.7p includes also an oblique component for the hypercharge flux, given by 

-^obliq. = -^(^-^1 ^ + dzi A dZ2)QY ■ (A.l) 

We choose to integrate this flux as, 

AoUiq. = ^^(^1 ^ + ^1 A dz2 - Z2 A dzi - Z2 A dzi)QY ■ (A.2) 

Wavefunctions can be solved by following the general procedure described in subsection 12.31 
We find that oblique fluxes do not affect the physics of the wavefunction in a qualitatively 
important way. 



A.l Wavefunctions for the 5m matter curve 

In this case the gauge covariant derivatives appearing in eq. ()2.25p get modified to 
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where we have defined ai = 2Mi — iq^a. The corresponding matrix B reads 
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with dimensionless eigenvalues given by the roots of the following depressed cubic equation, 
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In terms of these quantities the solution to eqs. ()2.33p for Mi > leads to the wavefunctions 
for the ground state of each of the four towers of fields localised in the 5a/ matter curve in 
the presence of oblique hypercharge flux. The expressions are given by eq. (j3.15p with, 
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A. 2 Wavefunctions for the 10m matter curve 



Wavefunctions for the 10m matter curve can be obtained from the 5m wavefunctions by 
applying the transformations (j3.25p . We obtain that the ground state wavefunctions for 
M2 > are given by the general expression eq. (I3.15|) with 



10m 

Pi 



1 



Re(a2) , 



i A 



10a 



+ Re(a2) 



2 a;^"^^ - Re(a2) 



Im(a2) , 



10j\/ 

P2 



10m 

Pa 



-Im(a2) , 



-Re(a2) 



(A.9) 
(A.IO) 



where we have defined 02 = 2M2 + igya- Eigenvalues are given by the roots of the cubic 
equation, 
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A. 3 Wavefunctions for the 5h Higgs curve 

In the presence of oblique hypercharge fiux the gauge covariant derivatives in the {u, f )-basis 
are deformed to, 
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where ah = 2Mi2 — igyct. The matrix B now reads 
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The eigenvalues are given by the roots of 



v± v_ 



■Ke{ah) = 



p= 1,2,3 



(A.16) 



and by Aq^ = 0. The corresponding eigenvectors are, 



so 



/1\ 





Vo/ 



■'P 



I 







p = 1,2,3 



(A.17) 



where now, 





(A.18) 



With this information we can solve eqs. (12.33P and look for the wavefunctions of the ground 
state fields localised in the hu curve. For simplicity here we only solve for the case Myi > 0. 
However, note from eq. (|A.15P that 
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(A.19) 



is independent of the oblique components of the hypercharge flux, and therefore the latter 

do not affect the chirality of the zero mode and therefore to double-triplet splitting. We find 

that in the presence of oblique hypercharge flux the ground state wavefunctions for each of 
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the four towers of fields in the hu are given by eq. (13.380 with = ^(\y^^ > 
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A. 4 Wavefunction overlaps 

The wavefunction overlaps are calculated as in the main text and we simply display the 
results here. 
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To calculate the leading contribution triple couplings we rewrite the polynomial prefactors 
in the wavefunctions in the coordinates u and w and then expand them collecting factors of 
WW, uu, uw and wu. Then making use of the integral p.l7p we can write this as 
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Where is the Heaviside theta functions, 5k = ka + kf, — k^ — k^ and we have defined the 
quantities 

1 



Pi 

P2 
P3 



P4 



^1^ + 2 



^5m ^^5m _p5M _^5Mj ^ J'^IOm +^10m _ ^IOm _ ^IOm 
p5M ^p5M +p5M +^5m^ + (^^IOm +^10m + ^IOm + ^IOm 
_p5M +p5M _p5M +p5Mj ^ |^_plOM ^^IOm _plOM ^^IOm 
_p5M ^^5m +p5M _^5m\ ^ /'.plOM ^^IOm + ^IOm _ ^IOm 



and 



1 



1 



kllApl"-\<^H]-p\"kl« + 



lie; 



-ai 



{i,w) 

■Pl"k^"u) + 





f ah 


< 






vl 




( OLh 


* 




^vl 


vl 



(A.23) 



B Delocalisation and normalisation of wavefunctions 



The local form of the wavefunctions (j3.14p exhibits exponential localisation along all four 
real directions in S. Localisation onto a complex curve is induced by the Higgs vev, while 
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localisation within the curve is induced by the flux. However the general form of the wave- 
function ()3.13p need not be fully localised in this way. Indeed the arbitrary holomorphic 
prefactor can include an exponential so as to cancel the flux induced exponential localisation 
along one real direction of the matter curve. This situation is rather general for instance in 
toroidal compactifications [9]. It is clear that such delocalisation cannot occur over the full 
complex curve since a holomorphic function cannot cancel the non-holomorphic flux induced 
localisation. The local form of the wavefunctions (j3.14p are defined as the appropriate linear 
combinations of the, possibly delocalised, wavefunctions such that their expansion around 
the origin begins with increasing powers of z. This is a valid form to take for calculating 
the triple overlap since then the localisation of the other wavefunctions ensures that only 
this local form of the wavefunction contributes to the integral. However, for calculating the 
normalisation of the wavefunction the possible delocalisation can affect the results. Indeed 
the normalisation integral can diverge as an integral over and must be cutoff at the KK 
scale of S which for a homogenous manifold matches the local scale R^. 

Determining whether a wavefunction delocalises along one real direction in this way is of 
course a global question and not answerable in general from a local perspective. Therefore the 
normalisation integral of the wavefunctions carries an ambiguity. In this appendix we study 
this ambiguity quantitatively. Our aim is to quantify how much the normalisation integral can 
change when a wavefunction delocalises in this way. We do this by plotting the normalisation 
integral as evaluated in (j3.19p for the different generations for different values of the flux and 
of K. We also plot for the same values of flux and R the normalisation integral but with the 
exponential localisation factor along the curve dropped, thereby modeling the delocalisation 
effect. Of course dropping the exponential factor causes the wavefunction to delocalise along 
the full curve while we expect only a possible delocalisation along one real direction. But 
then this half-delocalistion can be expected to be somewhere close to the geometric average 
of the fully localised wavefunction and the wavefunctions with the exponential localisation 
along the curve dropped. Dropping the exponential localisation along the curve means that 
we perform the normalisation integral by explicit integration perpendicular to the curve while 
for the integral along the curve we simply take the homogenous answer of iZ| , similarly to the 
Higgs case in (j3.46p . We plot the resulting normalisation factors in figure [H As can be seen 
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Figure 9: Normalisation factors for localised and delocalised wavefunctions for flux values 
VI = 5/6, V2 = 5/4, Ml = 1.6M, M2 = 2M, 7 = 0, e = 1/10. The first plot fixes M = 1 and 
varies R while the second plot fixes R = 1 and varies M. 

the ratio of the delocalised wavefunction to the localised ones does not depend strongly on R. 
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The dependence on M however is stronger which is expected given that the amphtude of the 
flux measures the degree of locahsation along the curve. For large flux values the delocalised 
normalisation and localised ones, particularly for the lighter generations, differ significantly. 
However in the actual models we study the fluxes are fixed to be of order one, which means 
that the wavefunctions are already spread out to a significant fraction of and hence, as 
can be seen in the plot, the normalisations do not differ significantly. This means that we 
can be confident that the normalisation of the wavefunctions that we use are correct up to 
order one factors even if the wavefunctions delocalise. To illustrate this we plot in figure [ini 
the same couplings as in figure [5] except now we take the geometric average of the localised 
normalisation and the completely delocalised one, this being a good measure of delocalisation 
along one real direction. We see that the results do not differ significantly from figure [5l Note 




5 10 15 20 25 



R 

Figure 10: Plot showing the coupling as in figure [5] but with the normalisation of the matter 
curve wavefunctions calculated using the geometric average of the localised and delocalised 
normalisations. 

however that if taking the delocalised integral to be of order is not accurate up to order 
one factors then these results may be modified. This can occur if the manifold is strongly 
inhomogeneous. In that case a more global calculation would be required to obtain more 
accurate results. 
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